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20.  ABSTRACT  (Continued) 

sidelobe  canceller  (SLC)  configuration.  Complex  Gaussian  channel  fluctuations 
are  assumed  for  the  disturbed  ionosphere. 

Three  types  of  output  SNR  analyses  are  conducted:  non-fading,  fading 
with  array  time  constants  short  enough  to  adapt  to  fading,  and  fading  with 
array  time  constants  too  slow  to  adapt  to  fading.  Most  detailed  analytical 
results  are  carried  out  for  one  or  two  antenna  elements  and  slow  fading  time 
constants.  However,  some  more  general  results  on  processing  gain  are  obtained 
for  a  linear  array  of  M  SLC  antennas  and  slow  fading  time  constants  where 
processing  gain  is  defined  here  as  the  ratio  of  output  SNR  with  the  array 
functioning  to  achieve  optimum  performance  to  output  SNR  with  no  adaptive 
array.  Curves  and  simple  analytic  expressions  are  given  relating  processing 
gain  to  a  suitable  decorrelation  parameter.  The  latter,  in  turn,  is  related 
to  the  rms  angle  of  arrival  fluctuations. 

Error  rate  expressions  are  developed  for  coherent  and  incoherent  FSK,  PSK, 
and  DPSK  modems  when  used  in  conjunction  with  an  SLC  adaptive  array.  While 
general  analytical  results  are  obtained,  numerical  results  are  obtained  for 
the  simplest  case  of  a  single  SLC  antenna  and  a  single  jammer. 
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SECTION  1 


EFFECT  OF  A  SPACE/TIME  VARIANT  CHANNEL  ON  ADAPTIVE  ARRAY  PERFORMANCE: 

SPATIAL  DECORRELATION  EFFECTS 

This  document  is  a  topical  report  under  Contract  DNA001-82-C-0273 
covering  the  period  June  21,  1982  to  March  15,  1983.  The  subject  matter  con¬ 
cerns  the  effect  of  a  space/time  variant  channel  on  the  ability  of  an  adaptive 
array  to  null  jammers.  Attention  is  concentrated  on  the  effect  of  spatial 
decorrelation  between  antenna  elements  (the  elements  themselves  are  assumed 
small  enough  not  to  be  affected  by  spatial  decorrelation). 

While  other  adaptive  array  types  are  considered,  more  analysis  and 
evaluation  is  carried  out  for  the  sidelobe  canceller  (SLC)  configuration.  For 
detailed  performance  evaluation  a  Gaussian  spatial  correlation  function  and 
complex  Gaussian  statistics  are  assumed  for  channel  fluctuations.  These 
assumptions  are  appropriate  to  transionospheric  propagation  with  saturated 
scinti 1 1 ation  . 

1.1  CONTENTS 

The  subsequent  report  is  divided  into  seven  sections  and  one  appen¬ 
dix,  Section  2  presents  the  configurations  for  three  basic  types  of  adaptive 
arrays:  the  fully  adaptive  array,  the  multiple  beam  array,  and  the  sidelobe 
cancel  1 er . 

Section  3  presents  the  equations  defining  the  weights  for  optimum 
steady  state  performance  and  the  corresponding  mean  squared  error,  signal-to- 
noise  ratio,  and  antenna  pattern.  The  space/time  channel  averages  and  statis¬ 
tics  needed  to  calculate  performance  are  presented. 
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Section  4  presents  detailed  calculations  of  optimum  performance  for 
the  case  of  a  two-element  adaptive  array.  As  in  the  other  calculations  in 
this  report,  only  narrowband  signals  are  assumed  so  that  no  degradations  due 
to  frequency-selective  distortion  are  calculated.  This  was  done  to  place  em¬ 
phasis  on  the  impact  of  spatial  selectivity  (i.e.,  decorrelatio-')  or  direc- 
tion-of-arri val  fluctuations  upon  performance.  The  space/time  variant  channel 
model  without  frequency  selectivity  consists  of  a  collection  of  channels  from 
each  transmitter  to  each  antenna  element  with  each  such  channel  modeled  as  a 
complex  time-variant  multiplier.  Three  types  of  analyses  are  conducted:  non¬ 
fading,  fading  with  array  time  constants  short  enough  to  adapt  to  fading, 
fading  with  array  time  constants  too  slow  to  adapt  to  fading.  A  two-element 
half  wavelength-spaced  array  is  considered  and  response  to  a  single  jammer 
evaluated. 

Section  5  in  large  part  presents  detailed  calculations  for  the  SLC 
(sidelobe  canceller)  array  that  parallel  those  in  Section  4  for  the  fully 
adaptive  array.  In  this  case,  however,  a  linear  array  with  a  single  antenna 
and  two  auxiliary  antennas  is  considered  together  with  the  response  to  two 
jammers.  At  the  end  of  this  section  some  analytical  results  on  processing 
gain  against  a  single  jammer  at  large  J/N  ( jammer-to-noi se)  ratios  are  ob¬ 
tained  for  a  linear  array  of  M  SLC  auxiliary  antennas.  By  processing  gain  we 
mean  the  ratio  of  output  SNR  with  the  array  functioning  to  achieve  optimum 
performance  to  output  SNR  with  no  adaptive  array. 

Section  6  carries  out  further  analysis  of  SLC  processing  gain  for  a 
single  jammer,  a  linear  array  of  M  auxiliary  antennas,  and  array  time  con¬ 
stants  long  compared  to  the  fading  time  constants.  With  the  aid  of  an  appro¬ 
priate  spatial  correlation  function,  formulae  are  obtained  for  processing 
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gain.  Plots  are  presented  of  processing  gain  as  a  function  of  the  decorre¬ 
lation  between  fading  on  adjacent  antennas. 

Section  7  presents  the  basic  analysis  for  the  more  general  case  of 
a  three-dimensional  grid  of  SLC  auxiliary  antennas  and  multiple  jammers.  A 
computer  program  was  written  to  calculate  processing  gain  for  this  general 
case  but  it  has  not  been  exercised  to  any  extent  yet. 

Section  8  develops  expressions  for  the  error  rate  of  modems  used  in 
conjunction  with  an  SLC  (sidelobe  canceller)  adaptive  array.  It  is  assumed 
that  the  propagation  medium  produces  decorrelations  in  time,  space,  and  fre¬ 
quency,  on  the  transmitted  signals.  However,  it  is  assumed-  that  the  adaptive 
array  time  constants  are  too  slow  to  adapt  to  the  fading  and  that  the  signal 
bandwidth  is  narrow  enough  to  avoid  significant  frequency-selective  fading. 
Thus  the  performance  degradation  suffered  by  the  array  is  caused  by  spatial 
selectivity  and  nonadaptation  to  the  fading.  The  analysis  is  concerned  with 
the  error  rate  performance  of  a  digital  modem  connected  to  the  output  after 
the  SLC  has  attempted  to  cancel  the  jamming  signals.  As  in  the  case  of  the 
adaptive  array,  two  cases  may  be  considered:  the  modem  time  constants  may  or 
may  not  allow  adaptation  to  the  fading.  In  this  section  we  consider  the  case 
of  a  fading  adaptive  modem.  In  order  to  obtain  adequate  error  rate  perform¬ 
ance  it  is  necessary  to  employ  coding,  interleaving,  and  forward  error  cor¬ 
rection  techniques.  The  raw  error  rates  computed  here  may  be  used  to  evaluate 
the  improvement  offered  by  interleaved  hard-decision  error  correction  coding 
techniques.  Some  numerical  results  on  error  rate  are  obtained  showing  the 
impact  of  spatial  decorrelation. 
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1.2 


CONCLUSIONS 


In  the  case  of  a  fading  adaptive  array,  it  is  possible  in  theory  * 
for  the  antenna  pattern  nulls  to  become  sharper  with  increasing  jammers'  power 
until  the  jammers  are  completely  nulled  and  the  output  SNR  is  independent  of 
jammer  levels.  This  same  behavior  exists  for  the  non-fading  non-di spersi ve 
channel.  However,  there  are  degradations  over  the  non-fading  case  even  with  a 
fading  adaptive  array.  Aside  from  fading  per  se,  which  produces  large  modem 
performance  degradations  over  the  non-fading  case,  there  is  degradation  due  to 
the  fact  that  the  fluctuating  angle  of  arrival  of  the  signal  and  the  jammers 
may  cause  jammers  to  apparently  come  close  enough  to  the  signal  direction  to 
prevent  nulling  at  any  reasonable  jammer  level. 

In  the  case  of  the  fading  non-adaptive  array,  the  situation  is 
quite  different.  For  the  fully  adaptive  array  it  is  shown  that  as  the  fading 
becomes  severe  enough  so  that  the  mean  field  disappears,  the  array  "shuts 
down",  i.e.,  the  adaptive  weights  go  to  zero.  With  any  non-zero  mean  field 
the  output  SNR  at  large  jammer  levels  becomes  proportional  to  the  input  J/S 
ratio,  indicating  that  a  perfect  null  cannot  be  formed  at  large  jammer  powers 
as  in  the  fading  adaptive  case.  The  sidelobe  canceller  does  not  shut  down  as 
the  mean  field  goes  to  zero.  However,  the  output  SNR  is  proportional  to  the 
J/S  ratios  at  large  jammer  levels  also.  If  the  decorrelation  between  auxil¬ 
iary  antennas  is  sufficiently  small  and  if  sufficient  number  of  antennas  exist 
in  relation  to  the  number  of  jammers,  the  SLC  array  can  provide  adequate  proc¬ 
essing  gain  against  the  jammers.  The  reader  is  referred  to  Figures  6.2  and 
8.2  for  calculations  of  the  impact  of  spatial  decorrelation  on  processing  gain 

n - 

Practical  limitations  prevent  reaching  complete  jammer  cancellation.  We  do 
not  discuss  these  limitations  here  since  our  emphasis  is  on  determining  the 
impact  of  spatial  decorrelation. 


and  error  rate  for  the  case  of  a  single  jammer.  With  only  one  auxiliary 
antenna  the  impact  of  spatial  decorrelation  becomes  significant  with  small 
amounts  of  decorrelation.  Thus  from  Figure  8.2  it  may  be  determined  that  when 
the  correlation  coefficient  between  adjacent  auxiliary  antennas  drops  from 
1.0,  perfect  correlation,  to  0.9985,  the  PSK  modem  analyzed  suffers  approxi¬ 
mately  a  3  dB  degradation  in  SNR  performance.  However,  as  the  results  of 
Figure  6.2  indicate,  the  addition  of  another  auxiliary  antenna  will  cause  a 
very  large  improvement  in  processing  gain  and  remove  the  SNR  degradation  for 
the  same  amount  of  decorrelation. 

While  a  computer  program  has  been  developed  to  handle  quite  general 
three-dimensional  antenna  configurations  and  multiple  jammers,  there  has  not 
been  time  to  generate  results  for  this  report.  Also,  the  theory  developed  for 
error  rate  is  valid  for  multiple  jammers  and  antennas  but  insufficient  time 
was  available  to  calculate  results  for  this  report.  The  results  obtained 
should  be  extended  to  include  the  effect  of  time  and  frequency  selective  fad¬ 
ing.  When  the  bandwidths  become  so  large  that  system  operation  is  impossible 
due  to  frequency  selective  fading,  the  system  must  be  modified  to  include 
adaptive  equalize  s  at  the  outputs  of  antennas. 
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SECTION  11 


SYSTEMS  TO  BE  ANALYZED 

Figure  2-1  shows  a  simplified  block  diagram  of  a  fully  adaptive 
array  showing  M  elements  and  M  complex  weights  connected  to  the  antenna 
element  outputs  followed  by  a  combiner.  For  simplicity,  all  frequency 
conversion  and  filtering  operations  are  not  shown.  The  combiner  output  c(t) 
is  fed  to  the  modem  and  to  the  adaptation  circuitry.  A  reference  signal  r(t) 
is  extracted  from  the  modem  and/or  combiner  output.  An  error  signal  e(t)  is 
formed  by  subtracting  the  combiner  output  from  the  reference  signal 

e(t)  =  r(t)  -  c(t)  (2-1) 

The  error  signal  and  the  antenna  element  outputs  are  fed  to  a  signal  processor 
that  implements  an  adaptation  algorithm  for  the  complex  weights  p^,  ...  p^. 

The  transmitter  nodes  are  labeled  0,  1,  ...  L  with  the  signal  s(t) 
being  transmitted  from  node  0  and  the  jammers  j^(t),  ...  (t)  from  nodes  1, 

2,  ...  L,  respectively.  At  the  k'"^  antenna  element,  the  received  process  is 
represented  as  the  sum  of  three  terms 

w^(t)  =  S(^(t)  +  ij^{t)  +  n|^(t)  (2-2) 

a  signal  term  S|^(t),  a  jammer  term  i|^(t),  and  an  additive  thermal  noise  term 
n|^(t).  The  jammer  term  i  (^(t)  may  be  represented  as  the  sum  of  contributions 
from  the  L  jammers  to  the  antenna  element 

L 

ik(t)  =  Z  (2-3) 

m=l 

where  jamming  signal  received  at  the  k^^antenna  element  from  the 

th  . 

m  jammer. 
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Figure  2-1  Simplified  Diagram  of  Fully  Adaptive  Array 


*  1 


The  signal  term  S|^(t)  may  be  expressed  in  the  form 

S|^(t)  =  ys(t-c)  hgj^(t,t)  dc  (2-4) 

where  hQ|^(t,5)  is  the  time  variant  impulse  response  of  the  channel  from  node  0 
to  antenna  element  k.  In  general  we  represent  the  time-variant  impulse 
response  from  node  m  to  antenna  element  k  as  h^|^(t,5).  The  corresponding  time 
variant  transfer  function  is 

^mk(^>^)  =  (2-5) 

In  terms  of  the  time  variant  transfer  function 

s^(t)  =  /S(f)  HQ^(f.t)  df  (2-6) 

where  S(f)  is  the  spectrum  of  the  transmitted  signal. 

The  received  jamming  signal  ^  mk  has  the  two  equivalent  ex¬ 
pressions 


(2-7) 

Vk>'> 

=  yjjf) 

^mk^^’^)  df 

(2-8) 

where  J  ( f )  is  the 
m 

spectrum  of 

The  combiner  output  is  given  by 


The  system  of  Figure  2-1  uses  an  array  of  wide  beam  antenna  ele¬ 
ments  or  subarrays.  An  alternate  approach  is  to  form  multiple  orthogonal 
beams  by  a  linear  transformation  (e.g.,  Butler  matrix)  and  apply  the  adap¬ 
tivity  on  the  beam  port  outputs  as  shown  in  Figure  2-2.  Q  _<  M  beams  are 
formed  by  separate  linear  transformations  on  the  M  subarray  outputs.  The  Q 
beamport  outputs  are  fed  to  adaptive  complex  weights  and  the  adaptation 
circuitry  is  the  same  as  in  Figure  2-1. 

Another  system  configuration  of  interest  is  the  sidelobe  canceller 
or  power  nuller  shown  in  Figure  2-3.  Under  normal  conditions  the  desired 
signal  appears  only  in  the  main  beam  of  the  signal  antenna  and  the  jammer 
appears  in  the  sidelobes.  The  auxiliary  antennas  perceive  the  jamming  signals 
and  by  weighting  the  outputs  of  the  auxiliary  antennas  and  subtracting  from 
the  signal  antenna  output  attempts  to  cancel  the  jamming  signal.  The  system 
of  Figure  2-3  can  be  cast  in  the  framework  of  Figure  2-1  by  identifying  the 
signal  antenna  output  as  the  reference  signal.  With  ionospheric  scintillation 
present,  the  transmitted  signal  may  appear  in  the  sidelobes  and  auxiliary 
antenna  outputs  and  the  jammer  may  appear  in  the  signal  antenna.  For  this 
reason  we  have  shown  paths  from  each  transmitter  to  each  antenna. 

In  this  report  we  will  confine  our  attention  to  the  systems  of 
Figures  2-1  and  2-3. 


13 
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Simplified  Diagram  of  Adaptive  Beam  Configuration 


Figure  2-3  Simplified  Diagram  of  Sidelobe  Canceller  System 


SECTION  III 


OPTIMUM  PERFORMANCE 


In  the  performance  analyses  we  shall  have  occasion  to  assume  either 
complete  adaptation  to  the  fading  or  no  adaptation  to  the  fading.  In  the 
former  case  it  is  assumed  in  effect  that  the  time  constants  of  the  adaptation 
circuits  are  sufficiently  small  for  the  system  to  follow  the  instantaneous 
fluctuations  in  the  time  variant  transfer  functions  latter 
case  the  adaptive  circuits  average  over  the  fading.  (We  neglect  non-ideali¬ 
ties  in  the  adaptation  process  and  circuitry.)  In  the  former  case  the  weights 
{P|^(t);  k=l,2,...M}  are  continually  adjusted  to  minimize  the  mean  squared 
error  signal  le(t) 1 ’  averaged  over  the  transmitted  signal  and  jammer  statis¬ 
tics  and  the  thermal  noise  fluctuations.  The  resultant  mean  squared  error, 
weights,  antenna  pattern,  SNR's,  etc.  will  be  time  variable  even  with  non¬ 
moving  jammers  and  signal  sources  because  of  the  fluctuations  of  the  channels. 
In  the  latter  case  the  system  will  adapt  to  changes  in  the  channel  statistics 
only. 

It  is  well  known  (e.g.,  see  [3.1])  that  the  solution  for  the 
weights  which  yield  the  minimum  mean  squared  error  is  given  by 


where 


P 


R 


P 


T 


(3-1) 


(3-2) 


R 


T 


(3-3) 
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in  which  T  denotes  matrix  transpose  and  we  have  suppressed  the  time  variable 
to  simplify  notation.  It  should  be  understood  that  the  averages  assume  the 
channels  "frozen". 

The  value  of  the  minimum  mean  squared  error  is  given  by  [3.1] 
Tep"  =  i  (3-5) 

or,  alternatively,  from  (3-1) 

Tip  =  pp  -  R.^  P.  (3-6) 

where  the  symbol  H  denotes  a  combination  of  matrix  transpose  and  complex  con¬ 
jugate  operations  ("Hermitian"  transpose). 

To  formulate  an  expression  for  the  antenna  pattern,  it  is  necessary 
to  specify  the  geometry  of  the  antenna  subarray  locations.  For  simplicity  we 
assume  a  linear  array  with  half  wavelength  spacing  between  subarray  centers. 
In  such  a  case  the  antenna  gain  pattern  is  expressed  as 

G(0)  =  G^(«)  ^  gj’rdc-Dsin* 

k=l 

where  Gj(()i)  is  the  antenna  gain  pattern  for  a  single  subarray  and  ^  is  the 
angle  measured  from  boresight. 
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We  may  formulate  an  expression  for  the  output  SNR,  p.  This  dis¬ 
cussion  of  SNR  is  for  the  fully  adaptive  array  and  must  be  modified  for  the 
SLC  (sidelobe  canceller)  configuration.  See  Section  5. 


p 


M 

L  Pk 


k=l 


S  P 


M 

k=l 


Pk<'k*"k> 


P”  ( Ni  ♦  J )  P 


where  is  the  identity  matrix. 


N 


(3-8) 


(3-9) 


C  *C 

Si  S2 

•  •  • 

C  *C 

^1  M 

s  *s 

2  ^1 

1S2P 

•  •  • 

s  *s 
^2 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

V^I 

c  *c 
^2 

«  •  • 

^1*^2 

•  •  • 

1^21' 

•  •  • 

i  *i 

2  M 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

« 

^M*^2 

•  •  • 

(3-10) 


(3-11) 


and  it  has  been  assumed  that  the  thermal  noises  in  the  different  antenna 
elements  are  statistically  independent  and  of  the  same  level. 


Since 


C  =  S  +  J  +  NI 


(3-8)  may  be  expressed  in  the  alternate  form 


(3-12) 


C  P 


P^  (NI  +  J)  P 


-  1 


(3-13) 


and  from  (3-1), (3-6)  we  obtain  the  further  equivalent  expressions 


P 

**  P^  (NI  +  J)  P 


(3-14) 


trp  -  |eM 
P^  (NI  +  J)  P 


(3-15) 


We  now  consider  the  evaluation  of  the  averages  required  to  deter¬ 
mine  £,  |e I p,  and  G((|i),  In  the  evaluation  we  shall  assume  that  the  signal 
and  jamming  processes  are  statistically  stationary.  Thus  [3.2] 


^(f)  S(0'  =  P(f)  «(f-i) 


(3-16) 


(3-17) 


where  P(f),  P^(f)  are  the  power  spectra  of  the  transmitted  signal  and  the  m^*^ 
transmitted  jamming  signal  and  «(•)  is  the  unit  impulse  function.  We  shall 
also  assume  that  the  various  jamming  signals  are  statistically  independent  to 
simplify  notation. 
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The  received  signal  crosscorrelations  are  given  by 


s,^s*  =  //  S*Tf)S(£) 

Using  (3-16) 

=  /P(f)  HQ^*(f.t)  HQ^(f,t)  df  (3-20) 

The  received  jamming  signal  crosscorrelations  are  given  by 
^  ^mk  ^mq 

•  r//v(f)  J„(0  dfdt  (3-21) 

Using  (3-18) 

V;=  2  df  (3-22) 

We  see  that  in  general  it  is  necessary  to  specify  the  joint  sta¬ 
tistical  properties  of  the  various  transfer  functions  to  evaluate  the  statis¬ 
tical  properties  of  the  optimum  solution  to  the  adaptive  array  problem.  In 
the  case  of  strong  ionospheric  scintillations  the  channel  transfer  functions 
become  complex  Gaussian  processes.  Then  their  joint  statistical  properties 
can  be  determined  from  knowledge  of  the  correlation  functicuis 


'^nmkq^^’^’’^’^^  ’  {  k’.q  =  r,2’, ! . 


This  correlation  function  is  more  general  than  the  mutual  coherence  functions 
normally  considered,  in  that  we  are  asking  for  the  time,  frequency,  and  spa¬ 
tial  coherence  functions  for  received  signals  corresponding  to  transmission  of 
carriers  from  two  different  locations  in  space  rather  than  one. 
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SECTION  IV 


TWO-ELEMENT  EXAMPLE  FOR  FULLY  ADAPTIVE  ARRAY; 

NON-FREOJENCY  SELECTIVE  CASE 

In  this  section  we  apply  the  results  of  the  previous  section  to  the 
simplest  non-trivial  example,  the  case  of  two  omni-directional  antenna 
elements  separated  by  a  half-wavelength  as  shown  in  Figure  4-1.  Such  a  con¬ 
figuration  is  able  to  discriminate  against  a  single  jammer.  We  confine  our 
attention  in  this  section  to  signaling  elements  that  are  sufficiently  narrow- 
band  that  no  frequency  selective  distortion  is  caused  by  the  disturbed  trans- 
ionospheric  channel. 

4.1  PERFORMANCE  WITHOUT  CHANNEL  DISTURBANCE 

It  is  instructive  to  compare  the  performance  before  and  after  the 
onset  of  the  channel  disturbance.  In  the  absence  of  a  channel  disturbance,  it 
is  sufficient  to  model  the  quiet  or  non-di sturbed  channel  by  a  simple  phase 
shift  for  the  case  of  narrowband  systems,  as  assumed  here.  With  the  desired 
transmitted  signal  on  boresight  we  may  use 


Ho^(f,t) 

=  =  1 

(4-1) 

H^^(f,t) 

=  1 

(4-2) 

^  gjTrsine 

(4-3) 

notation 

we  have  neglected  path  losses. 

This  results  in  no  loss 

of  generality  since  our  answers  are  expressed  in  terms  of  recei ved  power 
level s . 
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With  the  degenerate  transfer  functions  (4-1)  through  (4-3) 


w^(t)  =  s(t) 

+  j(t)  +  nj(t) 

(4-4) 

w^  ( t )  =  s  ( t ) 

+  j(t)  +  n2(t) 

(4-5) 

and  from 

(3-4) 

S+J+N 

C  = 

S.je-J-rSine 

S+J+N 

where 

s  =  ■ 

hTW 

J  =  ' 

ijTtTF 

(4-6) 


(4-7) 

(4-8) 


The  reference  signal  r(t)  is  taken  equal  to  the  signal  s(t), 
r(t)  =  s(t)  (4-9) 

so  that  [see  (3-3)] 


R 


T 


(4-10) 


Then,  using  (3-1)  we  find  that  the  optimum  weights  are 


C 


S+J+N 

-1 

S 

S+jg-ji-sine 

S+J+N 

S 

(4-11) 
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Upon  carrying  out  the  algebra,  we  find  that 


P  = 


4) 


where 


J  gjirsinej  ^  j 

J  (j_  -jttsine)  ^  j 


1 


8  = 


4(^)(f)sin^(|^ine)  +  2(^  +  +  1 


(4-12) 


(4-13) 


The  minimum  mean  squared  error  is  obtained  by  using  (4-12)  and 
(4-10)  in  (3-6).  Carrying  out  the  operations  and  normalizing  the  error  to  the 
strength  of  the  desired  signal, 


S  “  6(2(fj)  +  1) 


For  a  large  jamming  signal  (J  >>  S)  and  e  not  near  zero 


(4-14) 


where 


(4-17) 


i|)  -  n  s  i  n  6 

Upon  carrying  out  the  calculations  we  find  that 


P  S  P  =  e‘ 


S(|)‘4 


2(^)  sin*(|-  sine)  +1 


For  the  output  jamming  and  noise  power  we  have 

P^(NJ_+J)P  = 

^  (l-ej*)+l 


Again,  carrying  out  the  operations 


P^(N1+J)P  =  6*  N(|)' 

(2^+1) 

2(2^sin2{^sine)  +1) 

which  yields  the  output  SNR 

P^  S  P 

_  2S 

^2^  sin2(|5in*e)  +  l'' 

^  £^(NJ_+J)P 

N 

\  / 

We  note  the  limiting  cases 

P  - 

0 

TT 

"  7 

»  ♦  j§In  f'"' 

9 

-  0 

P  -  ^  Sin2(^sin0) 

for 

~  sin»(^sine)  >>  1 

J+N  J  e'^'*' 

^(l-eJ'^)  +  l,^(l-e-J’^)  +  lj 

Je"-^*  J+N 

(4-18) 


(4-19) 

(4-20) 

(4-21) 

(4-22) 

(4-23) 

(4-24) 
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For  e  -►  ir/2  the  jammer  is  completely  nulled  by  the  same  weights  which  maximize 
the  desired  signal  level,  while  for  e  0  the  jammer  is  located  on  boresight 
(as  the  desired  transmitter)  and  no  jammer  suppression  is  possible. 

It  is  intuitively  clear  that  when  the  reference  signal  equals  the 
desired  received  signal  s(t),  an  appropriate  measure  of  SNR  is  provided  by  the 
ratio  of  the  signal  strength  S  of  the  reference  signal  to  the  mean  squared 
departure  of  the  output  signal  from  the  reference  signal,  i.e.,  the  ratio 
S/  lep.  From  (4-13)  and  (4-14)  we  find  that 

S  4(^)(|)sin2(^ine)  +  2(^  +  |)  +  1 

,  _ ^  '  (4-25) 

W  2(^H) 


Using  (4-21)  in  (4-24)  we  find  that 


p  +  1 


(4-26) 


The  expression  for  the  antenna  pattern  G{^}  is  given  by  (3-7). 
Using  (4-12)  in  (3-7)  with  G^(4i)  =  1, 


G(«] 


6^(|)' 


(^)^l_e-jiTsine 


)+l 


jTtSin^, 


(4-27) 
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As  we  have  mentioned,  when  e  =  n/2  the  janimer  is  completely  can¬ 
celled  and  the  antenna  pattern  is  set  for  maximum  gain  at  (ji  =  0.  Collecting 
terms  in  (4-27)  and  normalizing  to  G(0)  for  e  =  ii/2. 


G(i») 


+1)  cos( 


7  sin,^)  -  -rr 


cosulsine  -  j  sin^,) 


2—  + 


(4-28) 
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4.2  FLAT  FADING  WITH  SPATIAL  DECORRELATION:  ARRAY  ADAPTIVE  TO  FADING 

The  first  case  we  shall  examine  with  channel  distortion  is  that  of  suffi¬ 
ciently  narrowband  transmission  that  essentially  flat  fading  results.  However 
it  is  assumed  that  the  antennas  are  far  enough  apart  that  there  is  some  de- 
correlation  between  the  perceived  jammer  and  signal  channel  fluctuations  at 
the  two  antennas.  We  also  assume  in  this  section  that  the  time  constants  of 
the  adaptive  circuitry  are  small  enough  to  allow  adaptation  to  the  fading. 

With  the  above  assumptions  we  may  express  the  channel  transfer 
functions  in  the  degenerate  forms 


lgi(T,t) 

=  901 (t) 

(4-29) 

Ig2(fit) 

902 

(4-30) 

gii(t) 

(4-31) 

=  gi2(t) 

(4-32) 

in  which  9|^p(t)  becomes  a  complex  Gaussian  process  in  the  limit  of  strong 
scintillation. 

The  received  processes  are  now  given  by 

Wi(t)  =  gg^(t)s(t)  +  g^^(t)j(t)  +  n^(t)  (4-33) 

W2(t)  =  go2(t)s(t)  +  gQ2(t)j(t)  +  n^{t)  (4-34) 
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and  the  crosscorrelation  matrix  by  [see  (3-4)] 


1901^5  +  +  N  g*^gQ2S  + 

’019*025  ♦  '’02l'5  *  I9l2l’''  *  " 

where  we  have  temporarily  omitted  the  argument  (t)  from  the  complex  channel 
gains  to  simplify  notation. 

The  reference  signal-input  process  crosscorrelation  vector  takes 
the  form  [see  (3-3)  and  use  (4-9)], 

^  =  gJiS  ,  gg^S  (4-36) 

To  solve  for  the  adaptive  weights  we  require  C^'^.  This  inverse  is 

found  to  be 

-(9019525  +  yngiz^e"'^’'^''"®)  l9oil*S  +  +  N 

where  the  determinant 

A  =  (  1901^5  +  |g^2l'J  +  (1902^5  +  Ig^^l*  J  +  N)  - 

^01^02^  ^  ^  i  (4-38) 
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Upon  carrying  out  the  algebra  in  (4-36)  we  find 


A  =  SJ  Ig, 


01'  '^12 


"  (L'^OiI’  *  l902lG  ^  ^  *  '’12'‘D  ■')  * 


(4-39) 


We  compute  the  weight  vector  via  the  matrix  operation  (3-1) 


P  =  C  R  -  - 

—  —  —  A 


(^*01  1^12!^  ■  ^11^12^02  ^01 

^902  I9lll"  “  ^11^12901  ®  ^  ^02 


(4-40) 


The  mean  squared  error  is  obtained  from  (3-6) 


S  -  R  P  = 


J  (soi  ISizI*  ■  9uSi2902  ♦  9 

-I  Joi^'Soa^J  /  .  .  *  q 

|9o2  I9iiI  911912901  ®  ^ ^  9 


*  SN 
01  ^ 


+  g*  SN 


(4-41) 


The  result  is 


It  is  worth  pausing  to  examine  the  expression  (4-41).  Assuming  the 

usual  case  that  the  jammer  is  much  stronger  than  the  signal  and  the  signal 
much  stronger  than  the  noise,  i.e.. 


J  >>  S  >>  N 


(4-43; 


with  the  average  strengths  of  the  channel  gains  normalized  to  unity,  i.e.. 


ISoir 


(4-441 


then,  except  where  the  term 


A(9,t) 


1  -  e 


•j„sine/^ll\/^02\ 

\yi2A9oi/ 


is  very  much  smaller  than  unity,  we  have  the  approximation 


1912!^ 


S  iQinl^  /  A(e,t) 


(4-46) 


Let  us  first  examine  the  form  this  expression  takes  when  the  spatial 
decorrelation  vanishes,  i.e.,  when 


Then  the  first  factor  becomes  the  numeric  2  and 


(4-47) 


(4-48) 


A(0,t)  =  4  sinM^  sine) 


(4-49) 


yielding 


-  SI  na 


-5 


We  see  that  the  dependenc-:  of  the  mean  squared  error  on  the  jammer 
angular  position  e  is  the  same  as  the  non-fading  channel,  but  the  level  of  the 
mean  squared  error  is  strongly  influenced  by  the  instantaneous  strength  of  the 
signal  channel,  i.e.,  of  the  received  signal  level.  In  a  deep  signal  fade  the 
mean  squared  error  will  become  much  larger,  as  one  would  expect.  However  this 
is  the  normal  impact  of  a  fading  signal  and  must  be  counteracted  by  modem 
design.  At  any  given  signal  level,  however,  the  adaptive  array  is  doing  its 
job  of  suppressing  the  jammer  via  the  third  factor  in  ^4-50). 


Returning  to  (4-461  and  assuming  that  there  is  a  small  amount  of 
spoLial  decorre ; ation ,  the  first  factor  in  (4-46)  will  reside  close  ti  the 
numeric  twu  on  the  average.  The  second  factor  is  the  same  as  in  the  case  of 
zero  spatial  decorrelation.  Thus  it  is  primarily  the  third  factor  which 
brings  in  the  effect  of  jammer  angular  position  on  mean  squared  error  and  we 
must  examine  this  term  to  assess  the  impact  of  spatial  decorrelation  on  the 
dependence  of  jammer  nulling  on  jammer  angular  position. 


We  consider  now  the  calculation  of  output  SNR  (3-8).  The  relevant 
matrices  are 


S  =  S 


^01^0? 


^01^02 


19(321  * 


[4-5ii 


NI  +  J  = 


j  +  N 


^11^1 2"^  e 


jTiSi  ne 


giig^^J  J.N 


(4-5?) 
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After  laborious  calculation  of 


S  P 

P  ~  u  u 

P”(NJ_+J)p" 

it  was  found,  again  that 


!e|^ 


(4-53) 


(4-54) 


Apparently  the  relationship  (4-54)  is  valid  for  the  non-frequency  selective 
narrowband  channel  case.  It  is  not  known  under  what  more  general  conditions 
(4-54)  may  be  val id . 


Using  (4-46)  we  find  that 


p  - 


5  (^)  A(6,t)  Pjg^^P  1^021" 


(h) 


<F'  *  '9l2l')  * 


(4-55) 


In  the  absence  of  jamming 


(|)  (l9oiM  ♦  I902l’) 


(4-56) 


Thus  the  adaptive  array  operation  produces  dual  diversity  with  square  law 
combining  in  the  absence  of  Jamming. 


For  strong  jamming 


^  A(e.t)  >>  1  (4-57) 


Using  Eq.(3-7)  and  the  adaptive  weights  (4-40)  we  find  that  the 


antenna  gain  pattern 


G(0) 


.g-JiT(sin,^-sin6)^*^^ 


-j  ttSI  n(t) 


01  ^02 


(4-58) 


The  ratio  of  the  gain  towards  the  signal  location  ^=0  multiplied  by  the  signal 
power  to  the  gain  towards  the  jammer  location  i))=6  multiplied  by  the  jammer 
power  is  sometimes  of  interest.  This  is  just 


SG(0) 

JOTeT 


S 

N 


N  (^12'®  ^  ^^ll)  (^Ul‘^^02) 

W  (^12~^ll)  ^01^12^^®’^^  ^^0l‘^^02  ^  ) 


(4-59) 


Assuming  no  spatial  decorrelation  but  still  scintillation 


^12  "  ^11 


(4-60) 


^01  ^02 


(4-61) 


and  (4-59)  simplifies  to 


G(0)  ^  S 

gITT  ■  J 


|g^^  1 2  sin  ^  (1,  sin  p  +  1 


cos  (  w  sin  |-) 


(4-62) 


which  differs  from  the  non-fading  ideal  channel  case  only  in  that  the  factor 
|g^^|2  accounting  for  the  jammer  power  fluctuation  is  present. 


4.3  FLAT  FADING  WITH  SPATIAL  DECORRELATION:  ARRAY  NON-ADAPTIVE  TO  FADING 


In  the  previous  section  we  studied  the  performance  of  the  two- 
element  adaptive  array  assuming  that  the  adaptation  circuits  can  adapt  to  the 
assumed  flat  fading.  Here  we  consider  the  other  extreme:  adaptation  circuit 
time  constant  larger  than  the  fading  time  constant.  In  such  a  case  the  cir¬ 
cuits  average  over  both  the  input  processes  and  the  fading  fluctuations.  The 
crosscorrelation  matrix  C  (4-35)  and  crosscorrelation  vector  R^  (4-36)  must 
then  be  changed  to 


C  = 


jnsme 


ISqiI^S  +  +  N  ^01^02^  ^  ^11^12'^^ 

^01^02^  ^  '''  ^ 


’1P12 


02' 


*12' 


(4-63) 


where  we  have  assumed  r(t)  =  s(t). 


(4-64) 


We  note  that  if  the  fading  is  sufficiently  severe  that 

9oi  “  9g2  ”  0  (4-65) 

as  in  the  case  of  complex  Gaussian  fading,  R^  will  vanish  and  the  optimum 
weights  will  go  to  zero,  shutting  off  the  array.  This  phenomenon  does  not 
occur  in  the  case  of  a  sidelobe  canceller  as  will  be  discussed  in  the  follow¬ 
ing  section. 


We  normalize  the  channel  strengths  and  define  correlation  coeffi 


cients  and  mean  as  follows 


I90ll^  =  =  1911^  =  1912^  =  1 


=  ^01% 

‘’0  '  “ZZIZZZZZIZI 

y  l9iil^  |9o2!" 


9*i9i2 

y  1^021  ^ 

'"o  ""  W  ""  ^ 

With  these  definitions 

C  +  1  oji'Sine 

oqS  +  p^J  e'J 

S  +  J  +  N 


S  +  J  +  N 

»5S  t  e'J'"'"* 


(4-66) 


(4-67) 


(4-68) 


(4-69) 


(4-70) 


R 


T 


(4-71) 


The  inverse  correlation  matrix  is  given  by 


--1  _  1 


S  +  J  +  N 


(pqS  +  DnJ  e 


jiisine\  1 

r  "  / 


-Lp 


+  o*Je 


-jirsinei 


S  +  J  +  N 


(4-72) 


where 


A  =  S^(l  -  IpqI’)  +  J^(l  -  Ip^P)  +  2SJ(1  -  Refpgpje-^*}) 

+  2(S+J)N  +  N*  (4-73) 


Using  Eq.(3-1)  the  optimum  weights  are  given  by 


mSN 


I  (I-Pq)  +  ^  (1-^1  +  1 


I  (l-p^)  +  ^  (l-p|  +  1 


:4-74) 


Note  that  when  there  is  no  spatial  decorrelation,  i.e., 

Oq  =  Pi  =  1  (4-75) 

the  optimum  weights  differ  from  the  non-fading  case  only  by  the  factor  m.  In 
fact  it  is  seen  that  when  (4-75)  applies,  the  channel  complex  gains  may  be  re¬ 
defined  as  part  of  the  corresponding  transmitted  processes  s(t)  and  j(t). 
Consequently,  antenna  patterns  and  SNR's  will  then  be  the  same  as  in  the  non¬ 
fading  case.  Mean  squared  errors  will  differ,  however,  because  of  the  reduc¬ 
tion  in  weights  by  the  factor  m. 

If  the  spatial  decorrelations  become  so  severe  that 

Pi  =  Pq  =  0  (4-76) 

then  (4-74)  shows  that  the  adaptive  weights  become  equal  in  value.  This  will 
produce  an  antenna  pattern  optimized  for  the  signal  in  the  absence  of  jamming. 
Thus  the  jammer  nulling  capability  has  been  completely  destroyed. 


Consider  another  special  case.  Suppose  that  there  is  no  spatial 
decorrelation  for  the  signal,  that  the  magnitude  of  the  jammer  spatial  decor¬ 


relation  is  close  to  unity  but  there  is  a  non-zero  angle,  i.e.. 


‘’0  = 
loi! 


1 

y 


(4-77) 

(4-78) 

(4-79) 


then  we  see  that  the  jamming  is  evident  (apart  from  the  factor  m)  only  in  a 
repl acement 

-IT  sin  e  Y  -IT  sin  e  (4-80) 

i.e.,  from  the  point  of  view  of  the  adaptive  array  the  angular  position  of  the 
jammer  has  been  changed  to  a  new  angle  i|/,  where 

n  sin  Tj,  =  ir  sin  e  -  y  (4-81) 

That  is  to  say  the  null  will  be  steered  in  error  by  the  angle  iji-e. 


The  mean  squared  error  is  given  by  Eq.(3-6) 


,  2  = 


„  |m|^  S^N 

A 


(1  -  Reiog})  +  ^  (1  - 


N 


Re{pieJ’’^’"®})  +  1 


(4-82) 
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Note  that  unlike  the  situation  wherein  the  array  adapts  to  the  fading,  the 
mean  squared  error  does  not  vanish  when  the  thermal  noise  vanishes: 


(1  -  Re{p„})  +  4  (1  -  Re{p,e'^’’^’"®}) 

_ ^ _ r _ z _ 

(1  -  IpqI^)  +  (j)\l  -  IP^r)  +  2(|)(1  -  Re{p*p^eJ’'^'"®}) 

(4-83) 


For  large  J/S  ratios  this  expression  further  simplifies  to 


-  l-21mp 


S 

J 


1  -  Refp^eJ’^^i^S 

IPiM 


(4-84) 


which  shows  the  mean  squared  error  increasing  as  J/S  increases.  This  may  be 
contrasted  to  the  case  of  fading  adaptation  where  the  error  reaches  a  minimum 
at  high  J/S  ratios.  The  problem  here  appears  due  to  the  term  ( J/S) *( 1- | p^  1  * ) 
which  vanishes  when  spatial  decorrelation  is  absent.  It  will  be  shown  that  the 
term  (l-|pj^|2)  is  directly  proportional  to  variance  of  the  angle  of  arrival 
fluctuations  due  to  the  ionospheric  scintillation. 
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SECTION  V 


SIDELOBE  CANCELLER  EXAMPLES 

In  this  section  we  apply  the  results  of  Section  3  to  the  case  of  a 
two-element  SLC  (sidelobe  canceller)  as  shown  in  Figure  5-1.  We  assume  suffi¬ 
ciently  narrowband  transmission  that  frequency  selective  distortion  may  be 
ignored.  As  in  the  case  of  the  fully  adaptive  array  of  Section  4  we  consider 
three  cases  in  order;  non-fading,  fading  with  array  adaptation  to  the  fading, 
and  fading  with  array  non-adaptation  to  fading. 

To  simplify  the  analysis  we  shall  confine  our  attention  to  the  case 
wherein  the  signal  level  received  in  the  sidelobe  canceller  antennas  is  small 
enough  relative  to  the  jammer  level  that  it  may  be  neglected. 

5.1  NON-FADING  CASE 

Using  the  above  assumptions  and  unity  gain  channels,  the  signals 
received  in  the  three  antenna  beams  of  Figure  5-1  are  expressible  in  the  form 


WgCt) 


Wj(t) 


Aq  (  9  )  j  ( t )  +  Ag(92)j2(t)  +  nQ(t)  (5-1) 

d  .  d  . 

-j2ir-sine,  -j2Tr-sin9„ 

A^{0i)e  ^  ^  jl(t)  +  Aj(e2)e  ^  j^lt)  +  n^{t)  (5-2) 


w^lt) 


A2(ej)e 


j2ir^ine^ 


j^(t) 


j2Tr-sinej, 

+  j2(t)  +  n2(t) 


(5-3) 


where  A|^(e)  is  the  complex  antenna  pattern  for  the  antenna  as  a  function 
of  the  angle  measured  relative  to  boresight  (see  Figure  5-1)  and  e^,  and 
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$2  are  the  angular  positions  of  the  signal,  jammer  1  and  jammer  2  transmit¬ 
ters.  The  distance  between  the  electrical  centers  of  each  of  the  sidelobe 
antenna  canceller  beams  and  the  main  antenna  is  d.  To  simplify  notation  we 
shall  assume  identical  canceller  beam  patterns 

A^(*)  =  A2(«)  =  A{*)  (5-4) 

As  discussed  in  Section  2,  in  the  case  of  the  power  nulling  system 
the  reference  signal  r(t)  is  set  equal  to  the  main  beam  output,  i.e.. 


r(t)  =  Wg(t) 


(5-5) 


Thus  the  crosscorrelation  vector  R  [see  Eq.(3-3)]  is  given  by 


=  I  Ul  Ul*  W  Ul* 

l_  ^0  1  ’  2  _ 


(5-6) 


Computing  the  averages  we  find  that 


jti/2 

j’l'p/2 

-jilii  /2 

"oi«  -'l  * 

-j'l'o/2 

^02® 

where  we  have 

used 

the  definitions 

On 

= 

Ao(en)  A*(e,) 

'll 

= 

4.  ^  sin  0^ 

and  J  i s  the 
n 

strength  of  j^(t)  • 

(5-7) 


(5-8) 

(5-9) 
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The  next  quantity  needed  is  the  correlation  matrix  C  [see  Eq.(3-4)] 


which  we  readily  find  to  be 


C  = 


jl'i  j't'o 


-Jill  -J'l'P 


(5-10) 


in  which  we  have  used  the  definition 


=r,  =  l''0<»n: 


(5-11) 


The  inverse  correlation  matrix  is  given  by 


-1  _  1 


+  G2J2  >  N 


j'l'i 

-G^J^e  -  G2J2e 


■jltli  "J't'o 

-G^J^e  -  G^J^e  G^J^  +  G^J^  +  N 


(5-12) 


where 


*  •  4  GiGjOiOjSin- 


(^) 


+  2N(G^Jj+  G2J2)  +  N*  (5-13) 


The  optimum  weights  are  obtained  from  Eq.(3-1)  as 


P  = 


j'l'i  j'1'2 

GfOi  +  G2J2  N  -(G^J^E  +  G2J2E  ) 


-J'l'i  "Jllio 

-(GjOje  '*  G^J^e  h  Gjjj  ♦  G^O,  *  N 


j'l'i/?  j'l'2/2 

^Ol'^l®  +  HQ2O2® 

-ji|i,/2  -j'l'2/2 

Hoidje  +  HQ2J2® 


(5-14) 
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Carrying  out  the  indicated  operations  we  find 


P 


j'i'i/2  j'l'o/2' 


~0'lip/2  -J4i 

^2^01®  ■  ^1^02® 


l/2\  / 

y 


j*./2  jtj,/2 

e  +  ^q2'^2^ 


+  HQ2J2® 


f 

) 


(5-15) 


The  minimutn  mean  squared  error  (minimum  output  power  in  the  present 
case)  is  given  by  Eq.(3-6).  After  carrying  out  the  operations  indicated  we 
find  that 


|eM  =  GqqS  +  N  +  2  N  D  (5-lb) 

in  which 


Gno  =  =  n=l,2  (5-18) 
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In  the  case  of  a  single  jammer  (5-16)  simplifies  to 


Gqq S  +  N  +  2N 


(5-19) 


The  first  term  in  (5-16)  is  the  signal  output  power  while  the  re¬ 
maining  terms  constitute  the  output  noise  and  jamming  level.  Thus  the  output 
SNR  is  given  by 


_  S  1 

P  fti 


00  N  1  +  2D 


(5-20) 


For  large  J/N  ratios  (iii^^  ^ 


IMA(e2)J  V|AQ(e2)  IMA(e^)  r-2Re{A*(0^)A(6^)AQ(e2)A*(e2)}cos 


Cxt) 


4|A(e)|‘|A(e2)rsin« 


(5-21) 


Note  that  D  is  directly  proportional  to  the  level  of  jamming  power  entering 
the  main  lobe  (i,.e.,  the  gain  of  the  sidelobes  of  the  main  beam).  If  we 
assume  the  ideal  condition  that  the  complex  gain  of  the  canceller  antenna  is 
set  equal  to  the  complex  gain  of  the  sidelobes  of  the  main  beams  and  these 
gains  are  independent  of  jammer  position 


Ao(ei)  =  Ay(02)  =  A(0^)  =  A(e2) 

Then  we  see  that  D  =  1  and 


p 


So' 

"TrT 


(A„ 


A) 


(5-23) 


Note  that  the  SNR  has  decreased  by  a  factor  of  3  over  an  unjammed  case. 


On  the  other  hand  if  the  sidelobes  of  the  main  beam  are  much  weaker  than 
the  gain  of  the  canceller  beams,  i.e., 

|Ao(en)|  «  A(0^)  :  n,m=l,2  (5-24) 

then  we  see  that 

P  -  So  I  ^1^0 1 

and  no  loss  in  output  SNR  will  occur  relative  to  the  unjammed  case. 

In  general,  at  large  J/N  ratios,  the  jammers  are  cancelled  but  the 
resulting  thermal  noise  level  will  vary  depending  upon  the  size  of  the  main 

beam  sidelobe  gains  relative  to  the  canceller  beam  gains. 
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5.2  FLAT  FADING  WITH  SPATIAL  DECORRELATION:  ARRAY  ADAPTIVE  TO  FADING 


In  this  section  we  assume  the  channels  from  the  two  jammers  and  the  sig¬ 
nal  to  the  three  elements,  nine  channels  in  all,  are  subjected  to  flat  fading. 
Moreover  we  assume  that  the  channels  from  a  given  transmitter  to  the  different 
antennas  may  fade  differently.  The  received  process  in  the  main  beam,  the 
reference  signal,  may  be  expressed  in  the  general  form 

Wo(t)  =  hQQs(t)  +  h^gj^(t)  +  h2Qj2(t)  +  nQ(t)  (5-26) 

where,  for  simplicity,  we  have  incorporated  the  effect  of  the  main  beam 
antenna  pattern  on  the  received  signals  within  the  definition  of  the  complex 
channel  gains  hQQ(t),  h^^(t),  and  h^^lt). 

Similarly  the  received  processes  in  the  two  sidelobe  canceller 

beams  are 


w^(t) 

+  h2^j2(t) 

+  nj(t) 

(5-27) 

W2(t) 

“  ^12 

+  h22j2(t) 

+  02! t ) 

(5-28) 

where  the  antenna  patterns  of  the  canceller  beams  and  any  phase  shifts  caused 
by  antenna  spacing  are  incorporated  in  the  :  m,n  =  1,2}. 
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Relating  the 

notation  in  the  present  and  previous  sections, 

hoo(^) 

=  Aq(0^3)  goQ(t) 

(5-29) 

hio(t) 

Ao(6i)  gjQ(t) 

(5-30) 

^20  ^  ^  ^ 

=  Aq(62)  92q(^) 

(5-31) 

h^j(t) 

-j4)i/2 

=  A^(ej)  e  9ii(t) 

(5-32) 

h2j(t) 

=  A^(02)  e  92i(t) 

(5-33) 

jifi  j^/2 

"  A2(0^)  e  9^2^^^ 

(5-34) 

h22(^) 

-  ®  922^^^ 

(5-35) 

where  represent  the  complex  channel  gains  that  have  been  normalized 

(without  loss  of  generality)  to  unit  average  squared  magnitude 

1^  ■  1  (5-36) 

Using  (3-2)  and  (3-3)  with  r(t)  =  WQ(t)  and  carrying  out  the 
indicated  averages 
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R 


^10^11^1  ^20^21'^2 


(5-37) 


^10^12'^!  ■*"  ^2o‘^22'^2 


C  = 


IhiiMJi  +  1^21^02  +  N 
'll"12''l  ■  ■■’21  *^22^2 


h*ihi2Ji  +  '121*^22^2 


+  h.„h*^J^  ^  I^22''^2  ^  ^ 


(5-38) 


where  we  have  assumed  that  the  two  jammers  are  statistically  independent  and 
of  strengths  8^,02. 

The  inverse  to  C  is  given  by 


,-l  1 


|h^2MJl  +  |h22l'J2 


-(h*jhj2'Ji  +  h*^h22J2) 


’^*^11^12'^!  ^  *^21^22^2  ^  ^  I*^21^*'^2 


(5-39) 


where 


a  =  I^inh22-h^2''21>'  '^iV  ( |  *+ |f’i2  * ^  1^21  > 1^22  I  *^^^2  ^ 


(5-40) 


The  optimum  weights  for  jammer  cancellation  are  then  [from  4-70)] 

P  = 


1 

1^2l'^l  ^  l^22l*'^2  ^  ^ 

-(^*1^12'^!  ^21^12'^2^ 

^10^11'^l'^^20^21"’2 

a 

^  ^21^22'^2  ^ 

Ihli^Jl  +  |h2il^J2  ^  ^ 

'^10'^12'^l‘^^20^22'^2 

(5-41) 
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Carrying  out  the  algebra 


(h*2h*i-  h^^h*2)(h^^h2Q-h2^h^Q)J^J2+  ^  10^12'^  1^'^'^2U^22'^2^ 

(5-42) 


Continuing  on  to  compute  the  minimum  mean  squared  error  (in  the 
present  context  this  would  be  the  minimum  array  output  power)  using  Eq.(3-6) 
we  find  after  considerable  algebra  that 


Vn  An  A'^IO'^^02^ 


.( 


|h 


10'  \H 


'20'  \  N 


leM  =  Ihgo^S  +  N  +  N 


(|hnl^+|hi?lMHNr/(|h2iM+|h22lOA 

(5-43) 


where 


^12  -  |N^^li22  -  ^12*^211' 


'01 


l'^10'^21  ■  ^20*^11 


4g2  =  1^10*^22  ■  ^20^12!* 


(5-44) 

(5-45) 

(5-46) 


The  first  term  in  (5-43)  represents  the  output  signal  power  while 
the  second  two  terms  represent  the  output  noise  and  jamming  power.  It  follows 
that  the  output  SNR  p  is  given  by 


(5-47) 
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A 


c 

■  i- 


which  becomes  for  large  J/N  ratios  ^ 


l^ol^ 


The  quantities  are  functions  of  time  due  to  their  dependence  on 
the  channel  fluctuations,  so  that  the  resultant  SNR  (5-48)  is  also  time  vari¬ 
ant.  An  interpretation  of  (5-48)  is  that  at  large  J/N  ratios  the  Jammers  are 
cancelled  at  the  output  but  the  resultant  thermal  noise  is  time  variable  and 
of  level 


(,  ^02  \ 
\  "^12  / 


(5-49) 


For  a  single  jammer  we  see  that 


/  +  Ihi2r)  +  1 


(5-50) 


which  at  large  J^/N  values  becomes 


I  *^00 '  ^  N 


(5-51) 


5.3  FLAT  FADING  WITH  SPATIAL-DECORRELATION:  ARRAY  NON-AOAPTIVE  TO  FADING 

In  this  section  we  will  assume  flat  fading  and  spatial  decorre¬ 
lation  as  in  Section  5.2  but  we  now  assume  that  the  adaptation  circuits  cannot 
follow  the  fading.  Section  5.3.1  considers  the  case  of  two  jammers  and  two 
SLC  antennas  while  Section  5.3.2  examines  one  jammer  with  M  SLC  antennas. 


5.3.1  Two  Jammers  and  Two  SLC  Antennas 

The  reference  crosscorrelation  vector  R  and  correlation  matrix  C 

now  become 


R 


^lO^ll'^l  ^  ^20^21'^2 

^10^12'^!  ^  ^20^22"’2 


(5-62) 


C  = 


‘^11^12'^!  ^  ^21^22'^2 


i  ^22  ^  *^2  *  ^ 


:5-53) 


Using  (5-29)  through  (5-36),  the  averages  in  (5-52)  and  (6-53)  may 
be  related  to  jammer  positions,  antenna  patterns,  and  fading  correlation 


coefficients  as  follows 

^10^11 

= 

AJ(ei) 

ji^/2 

^  pn,io 

(5-54) 

h^oh^l 

AJ(e2) 

jil'2/2 

®  P21,20 

(5-55) 

^10^12 

=  Ao(ei) 

A^(0j) 

®  ^12,10 

(5-56) 

^20^22 

=  AQ(e2) 

A*(82) 

®  '’22,20 

(5-57) 
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=  |Ai(6i)I^- 

(5-58) 

IS  I’ 

=  |Ai(e2)r 

(5-59) 

l^2l^ 

=  |A2(9i)M 

(5-60) 

l^'22l* 

=  |A2(e2)l* 

jl>T 

=  A*(ej)  A2(6^)  e 

Jlllo 

=  A*(e2)  A2(02)  e 

(5-61) 

h*ihi2 

^11, 12 

(5-62) 

^21^22 

P21,22 

(5-63) 

where  we  have  defined  the  correlation  coefficient 


^mn^pq 


mn  ,pq 


V 


|g  I  ^ 

'  pq' 


(5-64) 


The  correlation  matrix  inverse  is  given  by 


^  ^  *^^11^12'^!  ^  ^21^22'^2^ 
^  ^21^22'^2  ^  ^  1^21  ^*^2  ^  ^ 


(5-65) 
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where  h  is  the  determinant  of  C 


21^22!  ^ 


'*'  (1^22^*  1^22^^  ^  ^^12^*  ^^21^*  ~  ^21^22  *  ^11^12  ”  ^11^12  *  ^21^22  ^'^l'^2 


+  Ndh^^l^  +  |hi2l'’'^l  ^  1^21!'^  ^ 


(5-66) 


Using  Eqs.(5-54)  through  (5-63)  to  relate  the  averages  in  (5-66)  to 
the  antenna  gains  and  channel  correlation  coefficients 


a  =  G^2Si^^  ■  ‘^2l'"22^^  "  •'^21, 22!’^  ^2* 


[' 


+  1  G-«G,-+  G«- G,  -  2RG{Art»  ArtoA, -I  A'frtprti 


22  11  21  12  "''"‘"2r22'"lld2^21,22^12 


where 


^mn  =  l%nl’ 


A_„  =  A„(9m) 
mn  n'  m 


(5-67) 

(5-68) 

(5-69) 
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Expressions  for  the  optimum  weights,  mean  squared  error,  and  SNR 
may  be  obtained  as  in  the  previous  sections  but  these  do  not  simplify  as  in 

the  previous  case  and  we  do  not  present  them.  A  subsequent  memo  will  present 
results  based  on  computer  evaluation.  However,  to  obtain  some  understanding 
of  the  effect  of  spatial  decorrelation  without  using  a  computer,  it  is  suffi¬ 
cient  to  consider  special  cases.  Consider  first  the  case  of  a  single  jammer. 
Then  (5-52),  (5-53)  and  (5-65)  simplify  to 


R  = 


h^ohti  J 


(5-70) 


4  N 


hiih*2'  J 


^12^  J 


(5-71) 


,-l  1 


IhiiT  J  +  N 


■  ^11^12 


l^llT  J  +  N 


(5-72) 


where 


Using  (3-1)  and  (3-4),  the  optimum  weights,  minimum  output  power  and  output 
SNR  are 


(l'il2l*  *^10*^11  “  *^11^12  ^10^^  J  +  N 


{|hiil2  hjQh*2  -  '^10*^12^  ^  ^  ^lo'^12 


(5-74) 


(i)’ V  (t)"  V  (i) 


11 


lei*  =  IHoqI'S  +  N  +  N  ^  (5-75) 

i^y  V  (r)‘  '  *  ‘ 


.  sCs)'  “3*  (S 

p  -  IV ^  f3v  ■  .  rJN*  ,  ,  rJ> 


(S'  *  (S  “4  *  (S  “6  *  > 


(5-76) 


where 


Di  =  l^olM^l^  l^2l*  -  l^ol'  l^ll^2l'  • 


11 


1^10^12'  ^  ^  '^10^11  *  '^11^12  *  *^10*^12  ^  (^"^0 


^9  ~  1*^10^*  1^12^*  ”  ^  l*^inl*  1^11  I*  ■  I  (5-78) 


'10"12' 


'IQi  r'li 


10"11 


56 


2 


(5-79) 


^11 1  *  1^12  I  *  ~  I 


12  12' 


D4  -  D2  +  D3 


(5-80) 


and 


l^iol*  I^li2 
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1*^10 1  ^  |Aq(0j)M 


Examining  the  SNR  (5-63)  for  large  J/S  ratios  we  see  that 

-  S  ^3 

^  *  1*^001"  J  ‘  D7 


(5-8i; 


(5-82) 


(5-83) 


Thus,  as  in  the  fully  adaptive  array  (Section  4,3),  when  the  adaptive  array 
does  not  adapt  to  the  fading,  the  jammer  can  be  reduced  but  not  cancelled  and 
the  output  SNR  varies  as  the  input  S/J  ratio  for  large  J/N  ratios. 

The  "processing  gain"  against  the  jammer  (ratio  of  output  to  input 
S/J  ratio)  is  just 

0. 


g  = 


(5-84) 


1^00  I ^  D 


1 


since  the  input  signal/jamming  power  ratio  is 


(S/J), 


input 


l9ool  ^ 

Igiol* 


S 

J 


(5-85) 
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The  presence  of  spatial  decorrelation  causes  the  term  Dj  to  be  non-vanishing 
and  results  in  the  poor  asymptotic  behavior  (5-71).  To  see  this,  we  may  use 
Eqs.  (5-54)  through  (5-63)  and  (5-70)  to  express  the  coefficients  in  terms 
of  antenna  gains  and  channel  correlation  coefficients.  The  result  is 

*^1  ^  ^10^11^12  Q  ■  1^11,12!^  "  1^10, Ilf*  ■  I'^10,12l* 

+  2  Re  {  pfjjQ  Piu2  ^10,12  Q  (5-86) 

Expressions  follow  for  0^  and  (and  thus  D^)  which  also  exhibit 

the  vanishing  behavior  with  the  absence  of  spatial  decorrelation. 


*^2  ""  ^10*^12  [_^  ■  I‘’l2,10f*3  ^^11^10 

*^3  ^  ^11^12  Q  ■  1^11,12!*  3 


lPll,10l‘J  (5-87) 

(5-88) 


Thus  the  processing  gain  may  be  expressed  as 


1  - 


9  = 


f^ll.l2f 


1  - 


f'’ll,12l*  "  •‘’IQ,!! 


2  . 


•‘’10,12'*  ^  '^^^‘’11,10  ‘’11,12  ‘’10,12  ^ 

(5-89) 


An  even  simpler  case  worth  examining  is  the  case  of  a  single  can¬ 
celler  antenna  and  a  single  jammer.  In  this  case  we  assume  g^^  =  0,  =  0. 
Then  in  Eq.(5-75),  =  0^  =  0  and 


\ 


i 

t  I 
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(5-90) 


(k)  '^n  *  1 


"00  N  f0\2 


Again  we  observe  the  same  type  of  behavior  at  large  J/N  ratios 


5.3.2  One  Jaminer  and  M  SLC  Antennas 

We  now  derive  a  general  expression  involving  onl>  channel  correla¬ 
tion  coefficients  for  the  high  J/S  ratio  processing  gain  against  one  jammer 
when  M  s  delobe  cancelling  antennas  are  used.  With  a  single  jammer 

Wo(t)  =  hQQ(t)  s(t)  +  hjQ(t)  j(t)  +  n^lt)  (5-93) 

v^(t)  =  h^^(t)  j(t)  +  n^(t)  ;  m  =  1,2,...M  (5-94) 

The  reference  signal  crosscorrelation  vector  and  the  input  process  covariance 
matrix  are  given  by 


^10^11 


^10^12 


(5-95) 


C  = 


TFjYT*'  J  +  N 


hiihi2  J 


h|yi^  J 


Th^fTJ  +  N  ... 


^11 


Ih^  J+N 


(5-96) 
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cients 


In  terms  of  antenna  gain  patterns  and  channel  correlation  coeffi 


h 


h* 

10  Im 


ht  h, 
Im  Ip 


®  PmO 

A  (eJ  A*{eJ  e  p 


(5-97) 

(5-98) 


where  we  have  used  the  simpler  notation 

^Im.lO  "  ^mO 
^lm,lp  ”  ^mp 

and  defined  as  the  phase  shift  in  the  received  jamming  signal  at  the  m^^ 
antenna  relative  to  that  at  the  main  antenna  due  to  differential  path  lengths. 


(5-98) 

(5-100) 


Using  (5-98)  in  (5-96)  and  assuming  that  the  jamming  power  is  large 
enough  relative  to  the  noise  that  N  may  be  neglected  in  (5-96),  we  see  that 
the  typical  term  in  (5-86)  may  be  expressed  in  the  form 

c  =  (Apl.j)  A;(ei)  e''’'’'''’''"'’  .„p  ,  (5-101) 

where 

=  1  ;  m  (5-102) 
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As  a  consequence  of  (5-101)  we  may  factor  C  into  the  form 


C  =  J  r^^  ^ 


(N=0) 


(5-103) 


1^ 

where  ^  is  a  diagonal  matrix 


A«  = 


A*(epe 


A*(e^)e 


J  Y- 


A*(ei)e 


Jym 


(5-104) 


and  r  is  a  matrix  of  crosscorrelation  coefficients 


★ 

P12 


’IM 


Pl2 

1 


’l.M 


(5-105) 


Using  (5-97),  (5-95)  becomes 


R  =  J  Ag(9j) 


Jyi 

A^ (9j)6  ^10 

Jy^ 

^2^  ® 1 ^20 


Jyi^ 

PMO 


(5-106) 
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The  minimum  output  power  is  [Eq.(3-5)  with  the  reference  signal 

r(t)  =  WQ{t)] 

!e|’  =  I hgg  I  ’  S  +Th^F  J  +  N  -  C"^  R  c  (5-107) 


and  the  output  SNR 


1*^00  ^ 

p  = _ 

J  +  N  -  R^  c"^  R 


(5-108) 


Using  (5-85)  and  letting  N  0,  the  expression  for  the  processing  gain  at  high 
J/S  ratios  becomes 


1 


(S/J), 


input 


JT  A‘^  _r"^  (A^)‘^  R 


-  (5-109) 


Carrying  out  the  operations  indicated  in  (5-109)  we  find 


where 


1  H  ,-l 

1  -  p^Im  Pm 
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(5-110) 


(5-111) 
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SECTION  VI 


SOME  AJ  PROCESSING  GAIN  RESULTS  AGAINST  A  SINGLE  JAMMER 
FOR  THE  FADING  NON-ADAPTIVE  SIDELOBE  CANCELLER 

6.1  INTRODUCTION 

In  order  to  obtain  numerical  results  on  AJ  performance  using  the 
theory  developed  in  Section  5.3.2,  it  is  necessary  to  specify  an  SLC  antenna 
array  configuration  and  a  spatial  correlation  function. 

For  an  arbitrary  direction  of  arrival  and  an  arbitrary  antenna  con¬ 
figuration,  the  electrical  centers  of  the  SLC  antennas  will  not  all  fall  in  a 
plane.  However,  the  correlation  distance  along  the  direction  of  propagation 
is  very  much  larger  than  that  in  the  perpendicular  plane.  We  shall  assume 
that  for  the  size  antennas  of  interest  the  dimensions  along  the  propagation 
direction  will  be  small  enough  to  produce  negligible  spatial  decorrelation. 
As  a  result,  for  purposes  of  calculating  the  spatial  correlation  coefficients, 
it  is  sufficient  to  project  the  electrical  centers  of  the  antennas  into  a 
plane  perpendicular  to  the  direction  of  propagation.  The  distance  between 
electrical  centers  in  the  projected  antenna  are  used  in  the  computation  of 
correlation  coefficients. 

As  an  example,  consider  the  linear  array  shown  in  Figure  6-1.  The 
angle  of  arrival  of  the  jammer  (in  the  absence  of  scintillation)  is  e  relative 
to  boresight.  Antenna  numbers  1,2, ...,M  are  the  SLC  antennas  spaced  d  units 
apart  while  the  mainbeam  antenna  is  located  at  position  0.  The  projected 
antenna  is  shown  in  dashed  lines.  For  the  projected  antenna  the  elemehts  are 
spaced  d  cos  e  units  apart. 
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In  the  discussion  below  we  shall  initially  assume  a  general  spatial 
correlation  function  and  then  one  of  specific  form  appropriate  to  ionospheric 
scintillations.  Section  6.2  derives  expressions  for  the  AJ  processing  gain  of 
an  SLC  for  special  cases  of  interest  where  the  general  expressions  simplify  to 
a  greater  or  lesser  extent.  Section  6.3  presents  some  numerical  results. 
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1. 


6.2  DERIVATION  OF  AJ  PROCESSING  GAIN  FOR  SPECIAL  CASES 


According  to  Eq. (5-108),  the  output  SNR  for  a  single  jammer  is  given  by 


p 


I  ^00 1*  ^ 


Z-H—Z- - H - 1 

|h  J  +  N  -  R  C  R 


'10' 


(6-1) 


where  S  is  the  transmitted  signal  power,  J  is  the  transmitted  jammer  power  and 
N  is  the  additive  noise  power.  The  vector  ^  is  the  crosscorrelation  between 
the  main  beam  signal  and  each  of  the  SLC  antenna  output  signals,  and  £  is  the 
crosscorrelation  matrix  of  the  SLC  antenna  output  signals.  'In  the  absence  of 
the  SLC,  the  output  SNR  would  be 


1^0  s 

Po  ■  - - 

l^ol*  ^ 


(6-2) 


Defining  the  AJ  processing  gain  as  the  ratio  of  output  SNR  with  the 
adaptive  array  to  output  SNR  without  the  array,  we  find 

9„  ■  ^  ‘ _  (6-3) 

‘’o  1  -  R 

J  +  N 


Expressions  for  the  typical  terms  in  ^  and  £  are  [see  (5-93) 
through  (5-11)  modified  so  that  N  ^  0], 

R  =  J  AQ(e)  {A;(9)  e'^^"’  ;  m  =  1,2, ...M)  (6-4) 


C  = 


{Ap(e) 


-j(- 


*^m^ 


‘^mp 


+  6n,p  N  ;  m,p  =  1,2,...M}  (6-5) 
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where  A^{*)  is  the  antenna  pattern  for  the  antenna,  e  is  the  direction  of 
arrival  of  the  jammer  relative  to  boresight,  y  is  the  phase  delay  in  the  p^^ 
antenna  relative  to  the  main  beam  antenna  for  the  jamming  signal,  is  the 
crosscorrelation  coefficient  between  fluctuations  of  the  channels  from  the 
jammer  to  the  m^*^  and  p^^  antennas  and 
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mp 


1 
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m  =  p 

m  p 


(6-6) 


As  in  the  case  where  the  noise  N  was  assumed  zero  [Eq. (5-103)],  it 
is  possible  to  factor  C  as  follows 


C  =  J  A  r^  A 
—  —  — M  — 


(6-7) 


where  A  is  the  diagonal  matrix  [Eq. (5-104)] 


and  is  now 


A  =  {A^(e)  e  6|j,p  ;  m,p  =  1,2,...M} 


(6-8) 
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(6-9) 
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where 


mm 


IV»>I 


.th 


is  the  gain  pattern  of  the  m  antenna  evaluated  at  angle  e. 


(6-10) 


Using  (6-4)  and  (6-7)  we  find  that 


but 


where 


so  that 


^  R  =  -J-  A'^  r^^  (A^)"^  R 


rV'^  =  J  A*(p)  ^ 


M  10’  20’  ••• 


^m'2 


R^C-^  R  =  J  pj;  r;;^  P,^  |Ao(e)|* 


(6-11) 


(6-12) 


(6-13) 


Using  (6-13)  in  (6-3)  and  noting  that 


l^ol*  =  |Ao(e)M  =  G, 


0 


(6-14) 


we  obtain  the  general  expression  for  AJ  processing  gain  with  one  jammer, 

n  .  1 


’M 


1  - 


1  H  -1 


(6-15) 


1  +  N/JG, 
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We  now  consider  the  simpler  but  practically  interesting  case  in 
which  the  gains  of  the  SLC  antennas  are  all  equal.  Then  we  can  factor  r^^^  and 
P|^  in  the  form 


where 


(1  +  N/JGj^) 


(6-16) 


(1  +  N/JG^) 


(6-17: 


Im 


{£jjjp  »  i^»P  1»2»...  M  } 


1  pQp  *  m  —  M  } 


(6-18) 

(6-19) 


in  which 


and 


mp 


mm 


Pmp/(1  +  N/JGj)  ;  m  ^  p 

1  ;  m  =  p 

.  m—  1^2^...  M 


(6-20) 


(6-21) 


is  the  common  SLC  antenna  gain. 


Note  that  (6-18)  is  the  form  that  r^^  takes  when  the  receiver  noise 


is  set  equal  to  zero. 


becomes 


Using  (6-10),  (6-18)  and  (6-19)  in  (6-15),  the  AJ  processing  gain 


^-1  ^ 

^  _ 

(1  +  N/JG^)(1  +  N/JG^) 


(6-22) 


We  now  consider  the  further  simplifications  that  arise  when  the 
correlation  coefficient  p^p  takes  the  form 


=  p  =  R* 

^mp  m-p  p-m 


(6-23) 


This  form  arises  in  the  case  of  a  uniformly  spaced  linear  array  of  antennas  as 
shown  in  Figure  6-1.  For  this  case,  r^  and  take  the  forms 


where 


|_  ^1  ’  *^2  ’  *  •  •  3 


=  R^,  /(I  +  n/jg^: 


(6-24) 


(6-25) 


valid 


In  the  appendix  it  is  shown  that  when  the  forms  (6-24)  (6-25)  are 


''H  ;-i 


det  r 
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M+1 


det  r 


M 


so  that,  using  (6-26)  in  (6-22),  we  see  that 


(6-26 


(1  +  N/JG.)  (1  +  N/JGJ 


6.3  CALCULATIONS  FOR  A  LINEAR  ARRAY 

A  linear  uniformly  spaced  SLC  array  is  shown  in  Figure  6-1.  The 
antenna  locations  are  numbered  0,1,2,...M  corresponding  to  the  main  beam  and 
the  M  successive  SLC  antennas.  While  the  antenna  spacing  along  the  array  is  d 
units,  we  note  that  the  projected  antenna  spacing  is  d  cos  e.  Let  R(x)  denote 
the  correlation  coefficient  for  the  received  signal  fluctuations  (for  carrier 
transmission)  for  two  points  in  space  separated  by  x  units  and  lying  in  a 
plane  perpendicular  to  the  direction  of  propagation.  Then  we  readily  see  that 

Pmp  =  R(d[m-p]  cose)  =  R|^_p  (6-29) 

Assuming  the  case  of  saturated  ionospheric  scintillations,  Wittwer 
[6.1]  shows  that  the  spatial  correlation  function  R(*)  is  closely  given  by 


R(x)  =  exp  -  (6-30) 

*0 

where  Iq  is  a  correlation  "distance".  The  rms  value  of  the  angle  of  arrival 

fluctuations  a  may  be  related  to  ir.  by  the  following  equation  [Ref.  6.1] 
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n: 


so  that  (6-30)  may  be  expressed  in  the  alternate  form 


(6-31) 


R(x)  =  exp 
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Expressions  for  g|^  for  M  =  1,2,3  are  given  below 


1 

1 

- 

(6-36) 

1  -  0* 

1 

(6-37) 
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2a*  +  2a' 

'  -  2a*  (1  -  a*)(l  -  c*) 
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3a*  +  a' 

+  4a*  -  2a*  -  '2o*'*  -  2o‘*  +  4o"  + 

a“  -  30*'’  +  a"' 

ir 

a*)(l-a*) 

(6-38) 

It  appears,  but  has  not  been  proven  that 


9 


M 


_ 1 _ 

(1  -  aM(l  -  a')  ... 


(6-39) 


The  behavior  of  for  a  near  1  is  of  particular  interest  because 
it  describes  the  onset  of  degradation.  If  we  let 

a  =  1  -  €  (6-40) 

in  (6-39),  it  is  found  that 

-M 

9m  *  ^ 

2  M! 


From  (6-34)  we  see  that 


e 


(6-42) 


When  e  is  small  we  may  use  the  first  two  terms  in  a  Taylor  expansion  of  the 
exponential  in  (6-42) 


r~d  cose  * 
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It  follows  that  the  approximation  (6-41)  may  be  expressed  as 
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M  =  1,2.3 


(6-44) 


Figure  6-2  presents  plots  relating  g^ ,  g^  and  g^  in  dB  to  a  "decor 
relation"  parameter 

w  =  ^os  6  _  /T  ltd  cose 

— j-  -  - j -  Og  (6-45) 

using  the  exact  high  J/N  ratio  expressions  (6-36)  through  (6-38)  and  the 
Gaussian  shaped  spatial  correlation  function  (6-30).  The  dashed  lines  show 
the  small  decorrelation  approximation  (6-44).  Note  that  these  asymptotic 
approximations  which  are  straight  lines  on  a  log-log  plot  are  a  good  engineer¬ 
ing  approximation  for  values  of  the  decorrelation  parameter  y  as  large  as  0.3 
to  0.4  . 

The  improvement  in  performance  with  increasing  number  of  antennas 
may  be  understood  physically  by  noting  that  additional  antennas  allow  the  null 
to  be  broadened  and  thus  reduce  the  impact  of  the  angle  of  arrival  fluctua- 


DlCORRELATlUN  PAKAMFTtR  y 


Figure  6-2 

Decorrelation  Parameter  y  as  a  Function 
of  AJ  Processing  Gain  for  1,  2  and  3 
SLC  canceller  antennas  assuming  a  single 
jatimer  and  high  J/N  ratios.  Linear 
Antenna  Array  with  Uniformly  Spaced 
Antennas,  Gaussian  spatial  correlation 
functions  assumed. 

M  =  number  of  SLC  antennas 

Dashed  line  denotes  small  decorrelation 


where 

■?  =  direction  of  arrival  of  janmer 
d  =  spacing  between  antennas 
1  =  wavelength  of  carrier  frequency 
•.  =  rms  angle  of  arrival  "luctuation 

=  correlation  distance 

n  -WTj  -  -  - -jrj - *  jg-  -  -  - 

AJ  ‘ROCESSl’iG  GAri  :in  .JB) 
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It  is  of  interest  to  examine  the  anomalous  case  in  which  the 


spatial 

correlatiton 

function 

has  the  form 

R(x) 

=  exp 

U) 

(6-46) 

and  thus 

«m  = 

|m  1 

B '  ' 

(6-47) 

where 
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exp  /- 
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*^0  / 

(6-48) 

The  processing  gain  is  given  by  the  ratio 
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Note  that  this  result  is  independent  of  M  in  contrast  with  the  previous  exam¬ 
ple  and  in  apparent  disagreement  with  our  heuristic  arguments  on  the  benefits 
of  more  antennas.  However,  one  may  argue  that  the  reasoning  is  not  wrong  in 
this  case  but  ambiguous  because  the  correlation  function  (6-40)  corresponds  to 
the  non-physical  case  of  an  infinite  value  for  the  rms  angle  of  arrival  fluc¬ 
tuations.  Clearly  no  null  broadening  can  counteract  an  infinite  rms  angle  of 
arrival . 


SECTION  VII 


RESULTS  FOR  THE  THREE-DIMENSIONAL  SLC  ARRAY 
WITH  MULTIPLE  JAMMERS 

7.1  INTRODUCTION 

This  section  will  formally  state  the  equations  that  need  to  be 
programmed  to  compute  the  processing  gain  of  a  specified  three-dimensional  SLC 
array  when  subjected  to  L  jammers  arriving  from  specified  directions. 


7.2  FORMULATION 


The  antennas  are  assumed  arranged  in  a  uniformly  spaced  three- 
dimensional  grid  pattern  with  grids  parallel  to  the  x,  y,  and  z  axes.  Thus 
the  location  of  an  arbitrarily  selected  antenna  can  be  specified  by  a  radius 
vector 


Inipq  ■  i"<‘x  *  t'-l) 

where  i_,  are  unit  vectors  along  the  x,  y,  and  z  directions;  d  ,  d  ,  and 

X  y 

d^  are  the  uniform  spacing  between  antenna  elements  aligned  parallel  to  the  x, 
y,  and  z  directions,  respectively;  and  m,  p,  and  q  are  integers.  By  letting 
the  antenna  degenerate,  a  planar  and  simple  linear  array  can  be  modeled.  Thus 
the  formulation  includes  the  two-  and  one-dimensional  SLC  performances  as 
special  cases.  We  assume  that  m,  p,  and  q  range  over  the  values  (-M_,  M^) , 
(-P_,  P^)  and  (-Q.,Q^)  respectively. 


It  is  convenient  to  specify  an  antenna  and  its  output  by  the  sub¬ 
script  mpq.  We  assume  that  the  main  beam  antenna  is  located  at  the  origin  of 


where  j.(t)  is  the  jammer,  h,  (t)  is  the  complex  time  variant  gain 
l  j  »  j,,mpq 

characterizing  transmission  from  the  transmitter  location  to  the  mpq 

antenna  output,  and  n^pq{t)  is  the  thermal  noise  component  at  the  mpq  antenna 
output.  Note  that  we  have  assumed  no  signal  components  in  any  of  the  SLC 
antenna  outputs,  i.e.. 


0,mpq 


=  0 


mpq  ^  000 


(7-3) 


The  complex  transmission  gain  h  can  be  factored  into  the 

<,,mpq 

following  form 


h  (t)  =  A  (e  q  (t) 


(7-4) 


where  0re  the  spherical  angular  coordinates  characterizing  the  direc¬ 

tion  of  arrival  of  the  jammer,  A^pq(e»4')  is  the  complex  antenna  gain 
pattern  of  the  mpq  antenna,  ■rn,pq(9j^»i>i^)  ’s  the  phase  delay  of  the  signal 
arriving  at  the  mpq  antenna  from  the  jammer  relative  to  the  phase  delay  at 
the  000  position  t Ygoo^  ‘  mpq^^^  complex  time- 

variant  channel  gain  characterizing  the  path  from  the  jammer  to  the  mpq 
antenna. 


To  simplify  the  antenna  model,  a  default  option  should  allow  all  the  SLC 
antennas  to  have  the  same  pattern,  i.e., 

A^pq(6,(j.)  =  B(e,4i)  ;  mpq  ^  000  (7-5) 

while 

^000^®’^^  -  A{0,^)  (7-6) 
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To  simplify  subsequent  notation  we  assume  temporarily  that  a  one- 
one  mapping  has  been  established  between  the  triplet  values 

{mpq  ;  m  =  -1 ,0,1 , 

P  ~  -P_ , .  •  •  - 1 ,0 , 1 , .  •  •  , 

q  =  -1,0,1,...Q^  } 

and  the  integers  {n;n=0,l,  ...  K}  ,  the  triplet  corresponding  to  the  integer  n 
being  denoted  by  nipP^qf,.  We  assign  n=0  to  the  triplet  0,0,0  denoting  the  main 
beam. 

Thus,  rewriting  Eqs.(7-1)  through  (7-4), 

WQ(t)  =  s(t)  hgQ(t)  +  j^{t)  hjQ(t)  +  ...  +  j|^(t)  h|^Q(t)  +  nQ(t)  (7-7) 

w^(t)  =  j(t)  hj^(t)  +  j2(t)  h2^(t)  +  ...  +  j|_(t)  h|_r,(t)  +  n^(t)  (7-8) 

hgn  =  0  ;  n  t  0  (7-9) 

JYn^  ♦♦a) 

e  "  ‘  ‘  9,„(t)  (7-10) 


The  SNR  in  the  main  beam,  i.e.,  the  SNR  without  adaptation,  is 


%0  ^ 

‘^O  L 

E  +  N 

1=1 


(7-11) 


where  we  have  used  the  notation 


for  the  n  antenna  gain  in  the  direction  of  the  i  janmier  and  we  have  used 
the  normalization 


I  I*  =  1 


(7-13) 


The  SNR  after  subtraction  of  the  SLC  antenna  outputs,  i.e.,  the  SNR 
after  adaptation,  is  given  by 


^00  ^ 


.H  ^-1 


Z  J.  +  N  -  R"  C'"  R 


(7-14) 


so  that  the  processing  gain  is  given  by 


2 


E  J,  0^  *  H  -  r"  c"'  r 
1=1  *  *"  ~ 


where 


(7-15) 


—  "  ^  WgWjJ  ;  n  =  1,2,...,K  } 


(7-16) 


^  =  {  w^Wj;  ;  n,r  =  1,2,...,K  ) 


(7-17) 


Using  Eqs.(7-7)  through  (7-10)  we  can  express  ^  and  £  as  sums  of  matrices  each 
dependent  upon  a  single  jammer  (assuming  uncorrelated  jammers). 


P.nn  •  ~  1 ,2, . . ,  ,K  } 


-J  Crr(®, 

.♦.)  A;(.  ,♦.)  e  r  ,  .  n  .  . 


'■  "  ‘  p  +  N6 

inr  nr 

;  n,r  =  1,2 . K  }  (7-18) 


where 


1  ;  n  =  r 

0  ;  n  ^  r 

and 


(7-19) 


(7-20) 


is  the  spatial  correlation  coefficient  between  reception  at  the  n^  and  r^^ 
antennas  corresponding  to  transmission  from  the  position  of  the  jammer. 


Assuming  our  interest  is  in  SNR  or  processing  gain  we  can  normalize 

both  S  and  to  N,  i.e.,  set  N=1  in  the  above  expressions  and  regard  as 

the  J  /N  ratio, 
n 

It  remains  to  specify  the  correlation  coefficients  o.  „  and  the 
^  inr 

phase  delays  We  will  develop  expressions  for  these  terms  as  a 

1 1  Ht  X> 

function  of  the  antenna  element  positions,  the  arrival  angles  of  the  i*'" 
jammer  (e  ,♦  ),  and  a  specified  perpendicular  correlation  function  R|(x).  The 

X< 

latter  is  the  spatial  correlation  function  in  a  plane  perpendicular  to  the 
direction  of  propagation  of  a  received  signal. 


As  discussed  earlier,  in  the  computation  of  the  spatial  correlation 
coefficients,  it  is  sufficient  to  use  modified  antenna  element  positions  ob¬ 
tained  by  projecting  the  antenna  positions  on  a  plane  perpendicular  to  the 
direction  of  propagation  of  the  jammer.  Thus  if  d^^^  denotes  the  distance 
between  the  n*^^  and  r^"  projected  antennas  for  jammer  direction  of 

arrival , 

“<nr  = 
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Then  for  computation  of  it  remains  only  to  obtain  an  expression  for 

once  R|^(*)  is  specified. 

Consider  now  the  phase  delay  term  Assuming  that  the  jam¬ 
mer  is  in  the  far  zone  of  the  antenna  depends  only  on  modified 

antenna  positions  obtained  by  projecting  the  actual  antenna  positions  on  a 
line  from  the  jammer  to  the  main  beam  antenna.  Let  the  position  of  the 
main  beam  antenna  be  the  origin  of  coordinates  along  the  line  and  s^^^  be  the 
position  of  the  n'"^  projected  antenna,  then 

-  2"^  (7-22) 

where  x  is  the  wavelength. 

The  distances  d,  ,s„  may  be  obtained  from  the  antenna  positions 
tnr  an 

in  the  x,y,z  coordinate  system  by  means  of  a  rotation  of  coordinates.  Figure  1 

shows  the  relation  between  the  x,y,2  coordinate  system  and  the  x  ,y  ,z  coor- 

£ 

dinate  system.  The  z  axis  coincides  with  a  line  from  the  origin  of  x,y,z 

system  (the  main  beam  antenna  location)  to  the  jammer.  This  line  has 

spherical  angular  coordinates  The  y  axis  is  perpendicular  to  z  and 

£  £  £  £ 

in  the  plane  containing  the  z  and  z  axes.  Since  the  x  ,y  ,z  coordinate  sys- 

it  t 

tern  is  an  orthogonal  system  the  x  y  plane  is  perpendicular  to  the  z  axis.  It 

£  £ 

follows  from  our  construction  that  the  x  y  coordinates  of  the  antennae  are 

£ 

just  those  of  the  plane-projected  antennas  needed  to  compute  d^^^  while  the  z 

£  £  £ 

coordinate  provides  s^^^.  Specifically  if  {^n»yn’^n^  denote  the  coordinates  of 

t  h  ft.  i  i 

the  n  antenna  in  the  (x  ,y  ,z  )  coordinate  system. 
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tnr 


in 


V 


{xj  -  xj;  )  +  (yj  -  yf) 


(7-23) 

(7-24) 


The  relationship  between  (x^.y^.z^lj)  and  (x^.y^^.z^)  is  given  by 
"n"  =  "n  -  ^n 

cos<(,^  +  y^  sin(),^  -  (7-26) 

0 

z„  =  x„  sine  cos*  +  y„  sine  sin*  +  z„  cose  (7-27) 

n  n  i  •'n  n  "it 


From  (7-1)  we  note  that 


X 

y 


n 

n 


z 


n 


P  d 
^n  y 

Q  d 
^n  z 


(7-28) 

(7-29) 

(7-30) 


where  is  the  triplet  corresponding  to  the  n^*^  antenna  in  our  one-one 

mapping.  Using  (7-25)  through  (7-30)  in  (7-23)  and  (7-24)  we  see  that 


d*  „  =  d  (m  -m^)  sin*  -  d  (p  -p  )  cos*  * 
inr  |_  X  n  r  y'^^n  ^r 

[^^x^"'n""'r^  cos*^+  d^(p^-p^)  cose^  sin*^  -  d^(q^-q^) 

^n  =  "'n'^x  "  PnV’"®i  "’"♦i  "  Vz  "°^®i 


Sine  1’ 

(7-31) 


(7-32) 
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» 


As  discussed  earlier,  a  suitable  perpendicular  correlation  function 


is  given  by 

R^(x)  =  e  ^  (7-33) 

where  is  a  correlation  distance  parameter.  This  parameter  is  related  to 

a  ,  the  rms  value  of  angle  of  arrival  fluctions,  by 
0 


a 


6 


1 

ir 


*'0 


(7-34) 


The  above  formulae  allow  a  computation  of  processing  gain  for  a 

three-dimensional  array  subjected  to  noise  jamming  from  L  independent  sources 

arriving  at  angles  (e  ,41  )  ;  t  =  1,2, ...L. 

i  i 
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SECTION  VIII 


ERROR  RATES  WITH  A  FADING  NON-ADAPTIVE  SLC  ARRAY 

8.1  INTRODUCTION 

In  this  section  we  will  develop  expressions  for  the  error  rate  of 
modems  used  in  conjunction  with  an  SLC  (sidelobe  canceller)  adaptive  array.  It 
is  assumed  that  saturated  ionospheric  scintillations  nave  occurred,  causing 
the  propagation  medium  to  produce  complex  Gaussian  fluctuations  on  received 
signals.  It  is  also  assumed  that  the  adaptive  array  time  constants  are  too 
slow  to  adapt  to  the  fading  and  that  the  signal  bandwith  is  narrow  enough  to 
avoid  significant  frequency  selective  fading.  Thus  the  performance  degrada¬ 
tion  suffered  by  the  array  is  caused  by  spatial  selectivity  and  non-adaptation 
to  the  fading.  The  subsequent  analysis  will  be  concerned  with  the  error  rate 
performance  of  a  digital  modem  connected  to  the  output  after  the  SLC  has 
attempted  to  cancel  the  jamming  signals.  As  in  the  case  of  the  adaptive 
array,  two  cases  may  be  considered:  the  modem  time  constants  may  or  may  not 
allow  adaptation  to  the  fading.  In  this  section  we  consider  the  case  of  a 
fading  adaptive  modem.  In  order  for  adequate  error  rate  performance  to  result 
it  is  necessary  to  employ  coding,  interleaving  and  forward  error  correction 
techniques.  The  raw  error  rates  computed  here  may  be  used  to  evaluate  the 
improvement  offered  by  coding  techniques.  Section  8.7  presents  numerical 
results  on  error  for  the  simplest  case  -  one  SLC  antenna  and  one  jammer. 
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8.2  SYSTEM  DESCRIPTION 


The  system  to  be  analyzed  is  shown  in  Figure  8-1.  An  array  of  M 
auxiliary  antennas  with  adaptive  complex  weights  and  a  combiner  form  a  side- 
lobe  canceller  (SLC)  which  by  subtraction  attempts  to  cancel  the  L  jamming 
signals  appearing  in  the  sidelobes  of  the  signal  or  main  beam  antenna.  The 
"cancelled"  output  is  fed  both  to  a  digital  modem  for  extraction  of  the  output 
data  stream  and  to  a  device  which  implements  the  adaptation  algorithm.  This 
algorithm  attempts  to  set  the  complex  weights  p^,  p^.  Pj^  for  minimum 

output  power.  However,  the  adaptation  process  is  assumed  to  be  too  slow  to 
fol low  the  fading. 

The  combiner  output  c(t)  can  be  expressed  as 

L 

c{t)  =  s(t)  +  £  i^(t)  h^Q(t)  +  TiQ(t) 

-  E  pIe  *  n„(t))  (8-1) 

m=l  \ii=l  ' 

where  s(t)  is  the  transmitted  signal,  lij(t);  1=1 ,2 , . . .  ,L}  is  the  set  of  L 
jamming  signals,  is  the  thermal  noise  in  the  m^^  antenna  output,  and 
is  a  complex  time  variant  multiplier  characterizing  transmission  from  the 
jammer  to  the  m^*^  antenna  output.  We  may  rewrite  (8-1)  in  the  form 

L 

c(t)  =  s(t)  hf,(t)  +  Y.  h  (t)  +  n(t)  (8-2) 

^  ii=l  ^  ^ 
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where 


M 


h^(t)  =  h^o(t) 

-  S  p„  p„(t) 

m=l 

(8-3) 

n(t)  =  nQ(t) 

-  2  p„  p™(t) 
m=l 

(8-4) 

(8-5) 

Note  that  ij(t)  residual  jamming  component  after  cancellation  due 

to  the  i^^jammer,  while  n(t)  is  the  output  thermal  noise. 


The  transfer  functions  may  be  expressed  in  the  form 


V»t' 


«.)  e 


(8-6) 


where  the  angular  spherical  coordinates  of  the  jammer, 

Y  (6  is  the  phase  delay  of  the  received  jammer  signal  arriving  at  the 

m  i  i 

m^^  antenna  measured  relative  to  the  signal  antenna,  A^(e,^)  is  the  far  field 
pattern  of  the  m^^  antenna,  and  is  the  complex  time  variant  channel 

characterizing  the  path  from  the  jammer  to  the  m^^  antenna.  The  latter 
has  been  normalized  so  that 

■  1  <8-'> 

and  relative  phase  shifts  due  to  path  delay  difference  have  been  removed  and 
incorporated  in  the  phase  shifts  >  m=l,...,M}.  When  the  jammer  and 

signal  transmitters  are  far  enough  apart  and  saturated  fariing  occurs,  it  may 
be  assumed  that  the  path  gains  are  statistically  independent  complex  Gaussian 
variables. 
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<r, 


The  performance  of  coherent,  differentially  coherent,  and  in¬ 
coherent,  binary  and  quaternary  modems  will  be  evaluated.  In  all  cases  it  is 
assumed  that  the  error  rate  for  such  modems  with  a  non-fading  non-jammed 


channel  may  be  expressed  in  the  form 

=  ^(o)  (8-8) 

where  p  is  the  ratio  of  the  energy/bit  to  the  noise  power  density.  We  further 
assume  that  prior  to  detection  or  decision  operations  in  the  modem,  the  proc¬ 
essing  may  be  characterized  as  complex  linear  filtering.  The  jammers  are 
assumed  to  have  spectra  flat  over  the  bandwidth  of  this  linear  filtering  and 
to  have  filter  outputs  character! zable  as  independent  complex  Gaussian  proc¬ 
esses  when  channel  gains  are  fixed. 


As  a  consequence  of  the  above  assumptions  the  modem  error  rate 
conditioned  on  the  channel  complex  gains  is  given  by 


p 

‘  e 


(8-9) 


where  is  the  energy  bit  and  is  an  equivalent  noise  power  density 


Z  W  * 

i=l 


1  + 


M 

L 

m=l 


IPr. 


(8-lU) 


in  which  is  the  one-sided  power  spectral  density  of  the  real  jamming 
process  in  the  vicinity  of  the  sighal  carrier  frequency,  and  Nq  is  the  common 
one-sided  power  spectral  density  of  the  real  thermal  noises  at  the  antehna 
outputs.  may  be  expressed  in  the  form 


(8-11) 


Eb  '  jf/  "PW  '  If 

where  S  is  the  transmitted  power,  T  is  the  symbol  duration,  and  b  is  the 
number  of  bits/symbol . 

We  define 


'"ol’ 


s  iM-irrn- 

N  2b  '"o' 


(8-12) 


as  the  average  signal  energy/bit/thermal  noise  power  density  in  the  main  beam 
(W  is  the  receiver  bandwith  and  N  is  the  thermal  noise  power).,  and 


= 


(8-13) 


as  the  ratio  of  the  output  jammer  power  due  to  the  jammer  to  the  thermal 
noise  power  where  is  the  i'"  jammer  transmitter  power.  Equivalently,  n^^  is 
the  ratio  of  average  received  jammer  power  density  to  thermal  noise  for  the 
jammer.  With  these  definitions 


where 


Sol^O' 


Z  n  |z  1»  +  p 
£=1  ‘ 


(8-14) 


P  - 


1  ^  Z  iPn 


m=l 


(8-15) 
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8.3  ERROR  RATE  ANALYSIS  FOR  A  SINGLE  JAMMER: 

INDEPENDENTLY  FLUCTUATING  SIGNAL  AND  JAMMER  CHANNELS 

The  modem  error  rate  may  be  evaluated  by  averaging  the  conditional 

error  rate  (with  the  channel  gains  frozen)  over  the  channel  fluctuations, 

i  .e.. 


P  =  /  f(p)  W(p)  dp  (8-17) 

0 

where  W(p)  is  the  probability  density  function  of  p  in  Eq.(8-15).  We  consider 
now  the  evaluation  of  W(p). 

W(x)  may  be  obtained  as  the  derivative  of  the  probability  distribu¬ 
tion  function 

•W(x)  =  -gl  F(x)  (8-18) 

where 

F(x)  =  Pr  (p  <  x)  (8-19) 


Using  (8-15), 
Pr(p  <  x) 


Pr(So|Zor 


L 

1=1 


(8-20) 


We  consider  in  this  section  the  use  of  a  single  jammer  and  inde¬ 
pendent  fluctuations  for  the  signal  and  jammer  channels,  i.e.,  hjj  and  h^  are 
assumed  statistically  independent.  When  the  jammer  and  communicator  are 
separated  sufficiently,  the  statistical  independence  will  be  valid.  In  this 
case  we  need  tc  compute  the  canonic  probability 
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Pr(q  £  Q)  =  Pr(rQ  -  £  Q) 


(8-21) 


where 


k=0,l 


(8-22) 


and  ;  k=0,l}  are  independent  complex  Gaussian  variables.  The  probability 
density  function  of  r^^  is  given  by 


W.(r)  =  exp(  - -^)  U(s)  ;  k=0.1 


(8-23) 


where 


(8-24) 


U(r)  = 


r  >  0 


r  <  0 


(8-25) 


The  density  function  of  q  is  just  the  convolution  of  the  density 
function  of  r^^  with  that  of  -r^  or 

w(q)  =  /Wo(s)  W^(s-q)  ds 


Carrying  out  (8-26) 


aQ  +  a^ 


q 


(8-27) 


q 


The  density  function  of  ^  is  obtained  by  differentiating  (8-36) 
[see  (8-18)], 

Using  (8-17),  the  following  integral  yields  the  error  probability  for  the 
single  jammer  case: 

OB 

0 


An  alternate  expression,  which  is  sometimes  more  convenient,  may  be 
obtained  from  (8-38)  by  integrating  by  parts, 

CD 

%  =  /  r  ‘  -  T-rTTw'’”''!  Qn-Q'"  (6-W 

0 


The  following  represent  well-known  expressions  for  f(P), 


f(p)  =  ^exp  (-  |) 

;  incoherent  FSK 

(8-40) 

1  -P 

f(p)  =  -^e 

;  binary  and  quaternary  DPSK 

(8-41) 

f(P)  =  QC^) 

;  binary  and  quaternary  PSK 

(8-42) 

f(P)  =  Q(/— ) 

;  binary  coherent  FSK 

(8-43) 

where  FSK,  DPSK,  and  PSK  denote  binary  frequency  shift  keying,  differential 
phase  shift  keying,  and  phase  shift  keying,  respectively,  and 
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Q  )  =  /"  —  dY  (8-44) 

J 

/2p 

Noting  that 

1-^  exp  (- ^)  ;  incoherent  FSK 

j  exp  (-P)  ;  binary  and  quaternary  DPSK 

^  exp(-p)  ;  binary  and  quaternary  PSK 

exp(-  y)  ;  binary  coherent  FSK 

2/p  ✓2x 

and  using  (8-47)  in  (8-39),  we  find  that 

00 

'  iff  is C'  ■  ^ ^ 

0 

Using  the  change  of  variable  y  =  /2p, 

00 

Pg  =  y*  ^  exp(-  ^)  1^1  -  |:;yrp7g-  exp(-PyV2)  j  dy  (8-50) 
The  first  term  may  be  integrated,  yielding 


(8-45) 

(8-46) 

(8-47) 

(8-48) 


0 
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Upon  the  change  of  variable  y/R/2  z 


P 


e 


1 

2 


exp 


TTP" 


dz 


PSK 


More  generally,  since  (see  page  314  of  Ref.[8.1j) 


f  exp(-xxji 

Jq 


It  e^  Q(/^) 


we  see  that 


P 


e 


PSK 


Using  the  asymptotic  expansion 


Q{x) 


x/2^  l- 


we  deduce  the  approximations 


P 

e 


1  n  R  \ 


R  «  1+P 


R+P 

4 


R  «  1,  P  <<  1 
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(8-52) 


(8-53) 


(8-54) 


(8-55) 


(8-56) 


that 


In  an  exactly  analogous  procedure,  using  (8-48)  in  (8-39)  we  find 


P 


e 


H-2P 

2R 


coherent  FSK  (8-57) 


We  now  consider  binary  and  quaternary  DPSK  and  insert  (8-46)  in 
(8-39),  resulting  in  the  integral 


P 

e 


-p 

e 


1 

1+pR 


dP 


(8-58) 


This  integral  can  be  expressed  in  terms  of  the  exponential  integral 

function 


El(x)  = 


(8-59) 


as 


\ 

2 


e 


x4-P 

R 


1  +  P 
R 


DPSK  (8-60) 


With  the  3  dB  difference  for  incoherent  FSK 


1  - 


l-t-2P 
1  ^  2R 
7^  ® 


F 


1 

7 


incoherent  FSK 


(8-61) 
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8.4  ERROR  RATE  ANALYSIS  FOR  MULTIPLE  JAMMERS: 

INDEPENDENT  FLUCTUATING  SIGNAL  AND  JAMMER  CHANNELS 

In  this  section  we  consider  the  evaluation  of  error  rate  for  multi¬ 
ple  jammers.  Equations  (8-17)  through  (8-20)  are  applicable  to  the  multiple 
jammer  case.  Computation  of  (8-20)  requires  an  evaluation  of  a  probability  of 
the  generic  form. 


L 

H(Q)  =  Pr(r„  -  Y.  1  Q)  (8-64) 

1=1  ‘ 

where 

;  m  =  2,1,...  M  (8-65) 

and  Uj  is  a  complex  Gaussian  random  variable.  In  this  section  we  shall 
assume  that  the  jammers  and  communicator  are  mutually  separated  far  enough 
apart  that  the  set  {bj^(t)  ;  i=0,l,...L)  are  statistically  independent.  Then 
the  set  ;  i=0,l,...L}  may  be  taken  as  statistically  independent. 

We  shall  first  find  the  p.d.f.  (probability  density  function)  of 
L 

T  =  fQ  -  '"i  (8-66) 

1=1 

by  transforming  its  c.f.  (characteristic  function).  This  is,  if 

C(s)  =  e"^^  =  /w(q)  e’^^  dq  (8-67) 

is  the  c.f.  of  q  (we  have  generalized  to  a  complex  argument  s),  then  the 
p.d.f.  of  q  is  given  by 
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(8-68) 


f 


W(q)  =  J C(s)  ds 

Equations  (8-67)  and  (8-68)  are  two-sided  Laplace  transform  pairs.  The 
integral  in  (8-68)  is  a  contour  integral  closed  around  the  left  half  plane  for 
q  >  0  and  around  the  right  half  plane  for  q  <  0. 


Because  of  the  assumed  independence  of  the  set  the  c.f.  of  q 

factors  as  follows 


where 


C(s) 


L 

-s  (rg  -  r  r,) 

e  =  C(.(s)  n  c.(-s) 

^  1=1  ^ 


(8-69) 


t  =  0,1,...L  (8-70) 


is  the  c.f.  of  rj.  This  c.f.  is  given  by  [see  (8-23)  and  (8-24)] 


Ci(s) 


<D 


-sx  1  ,  X  >  . 

e  —  .exp(  -  —  )  dx 

1  1 


1/a, 

rriTi- 


(8-71) 
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Thus  using  (8-71)  in  (8-69) 


C(s) 


1/a, 


_ ^ 

s  +  1/a, 


0 


L 

n 

t=i 


( 


(8-72) 


We  desire  H(Q),  (8-64),  which  is  the  d.f.  (distribution  function) 
of  q  since  from  this  d.f.  the  d.f.  of  p.  F(x),  Eqs.(8-19)  and  (8-2U),  may  be 
trivially  determined  and  thereby  the  desired  error  rate  via  the  integral 


P 

e 


J  F(p)  (-f(p))  dp 
0 


(8-73) 


where  f(p)  is  given  by  Eqs.  (8-39)  through  (8-48). 

We  may  obtain  a  contour  integral  representation  of  H(Q)  as  follows 

^  Jr 

H(Q)  =  f  W(q)  dq  =  ^  y  ^  ds  (8-74) 

where  in  carrying  out  this  integration  the  path  of  integration  along  the  j 
axis  is  indented  to  the  right  of  the  pole  at  the  origin  and  the  contour  closes 
in  a  large  semi-circle  in  the  left  half  plane. 


Since  the  counter  includes  only  the  poles  at  s  =  -I/Sq  and  s  =  0, 
from  Cauchy's  residue  theorem. 


H(Q) 


s=-l/a 


or 


H(Q)  =  1 


■(!, «) 


0 


-Q/ar 


+  Res{^^  e^^l 


Q  >  0 


(8-75) 


s=0 


(8-76) 
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8.5  ERROR  RATE  ANALYSIb  FOR  MULTIPLE  JAMMERS: 

CORRELATED  JAMMER  CHANNELS  ONLY 

In  this  section  we  generalize  the  previous  analysis  to  the  case  in 
which  the  fading  on  the  various  jammer  channels  is  correlated  but  the  signal 
channel  is  still  uncorrelated  with  the  jammer  channels.  First,  we  wish  to 
compute  H(Q),  the  d.f.  of  the  random  variable  q  [see  Eqs.(8-64)  -  (8-66)]. 
This  d.f.  may  be  computed  from  (8-74)  once  the  c.f.  C(s)  is  known.  Under  the 
assumption  that  the  signal  channel  fluctuations  are  independent  of  the  jammer 
channels'  fluctuations,  C(s)  factors  into  the  product 


C(s) 


(8-84) 


The  first  average  has  already  been  computed  [see  Eq.(7.51)].  We 
consider  the  evaluation  of  the  second  average  which  is  the  c.f.  of  the  random 
variable 

L  L 

^  =  L  f'l  =  E  I**,!*  (8-85) 

1=1  ‘  t=l  ’■ 


Comparing  Eq.(8-20)  with  Eqs.(8-64)  and  (8-65)  we  can  make  the  identification 


u 


i  .e.. 


(8-86) 


L 

d  =  E  xn  Iz.l* 
1=1  ‘ 


(8-87) 
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The  random  variable  d  is  a  quadratic  form  in  correlated  complex 
zero  mean  Gaussian  variables.  Turin  [8.1]  has  evaluated  the  characteristic 
function  of  a  quadratic  form  in  complex  Gaussian  variables, 


d  =  /  qz 
as 


=  e  -  defi  X  +  sM  ^  J 


(8-88) 


(8-89) 


where  ^  is  the  matrix  defining  the  quadratic  form  and  M  is  the  moment  matrix 
of  the  random  variables.  In  our  case  ^  is  the  diagonal  matrix 


£  =  X 


where 


0  ...  0 

n  2  ...  0 

0 

0  n, 

and 

M  =  {  z*z^  ;  l,n  =  1,2,...L  1 

so  that 

"  detL  X  +  sM  N  J 


N  = 


(8-90) 


(8-91) 


(8-92) 


(8-93) 
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Using  (8-89)  and  (8-71)  in  (8-84) 


C(s) 


1/a, 


s  +  1/a, 


_ 1 _ 

det[  I  -  sx  M  N  ] 


(8-94) 


The  d.f.  H(Q)  is  thus  given  by  (see  (8-74))  the  contour  integral 

1  f  1  1 

~  JVJ  J  s  +  1/aQ  det[  j_  -  sx  M  N  ]  ?  (8-95) 

-j" 

It  will  be  shown  below  that  the  poles  of  the  second  factor  in  (8-74)  are  all 
in  the  right  half  plane.  Thus,  using  Cauchy's  residue  theorem. 


H(Q)  =  1-  - - e  ^  (8-96) 

det  C  I  +  ^  M  N  ] 

-  a^^  -- 

Using  the  definition  (8-32),  (8-34),  and  (8-81), 

where 


Rj  0  ...  0 

0  ...  0 

:  ■.  u 

0  0 


(8-98) 
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A  single  integral  expression  for  P  is  obtained  by  substituting  (8-97)  in 

^  • 

(8-73)  and  using  the  appropriate  expressions  for  (-f(p)).  Thus 


2  2  J 


e-(Pn)x 

det  LI  +  xM  RJ 


binary  and 
quaternary  OPSK 


(8-99) 


=  i  .  /■  J-  J- 
2  y  /77  /27 


-(P+l)x 


/jf?  det  [J,  +  xM  R] 


binary  and 
’quaternary  DPSK 


(8-100) 


Alternatively  we  note  that  the  denominator  of  (8-93)  is  an 
order  polynomial  in  sx  and  may  be  factored  to  explicitly  show  the  poles  of 
♦(s).  If  the  eigenvalues  of  the  matrix  product  M  are  denoted  by 
...LI  then 


det  [I  +  xM  N]  =  n  (1  -  sxx  ) 
“  “  ~  t=l 


(8-101) 


Due  to  the  positive  definite  nature  of  M  Ji,  all  the  eigenvalues  are  positive 
numbers  and  thus  the  poles  of  ♦(s)  are  in  the  right  half  plane.  Using  (8-101) 
in  (8-94)  we  arrive  at  the  following  expression  for  the  c.f,  of  q 


c(s)  = 


I,  ( 

V 


-  1/xx, 


(8-102) 


which  is  identical  to  the  form  (8-72)  in  the  case  of  independently  fluctuating 
jammer  channels. 
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It  follows  that  the  d.f.  H(Q)  is  given  by  an  analogous  equation  to 


(8-76) 


L  a„  -  Q/a. 

H(Q)  =  1  -  n  — -t~t-  e  ° 

t=l  ^  I 


(8-103) 


Then  using  (8-77)  and  (8-79),  the  c.f.  of  *  becomes 
F(x)  =  1 


(A 


-Px 


where 


St  = 


(8-104) 


(8-105) 


Finally,  the  error  rate  expressions  are  identical  to  (8-82)  and 
(8-83)  with  the  replacement  Sjj=Rj. 
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8.6  ERROR  RATE  ANALYSIS  FOR  MULTIPLE  JAMMERS: 

GENERAL  CORRELATED  CHANNEL  CASE 

We  consider  here  the  general  case  of  mutually  correlated  communi¬ 
cator  and  jammer  channels.  The  d.f.  H(Q)  can  be  computed  from  the  c.f.  C(s) 
via  (8-74).  We  may  use  the  results  of  Turin  [8.1]  again  to  compute  C(s), 
since  the  latter  is  the  c.f.  of  a  quadratic  form  in  complex  Gaussian  vari¬ 
ables,  i.e., 

i 


0 


0 


0 


Q  = 


-xn ,  0 


-xn^ 


-xn, 


(8-109) 


From  (8-74)  we  see  that 

H(Q)  =  ^  j  s  det  Lj_  +  sM  j  (8-110) 

left  half  plane 

where  the  pole  at  the  origin  is  included  in  the  left  half  plane  by  indenta¬ 
tion.  In  the  case  of  no  correlation  between  the  signal  and  jammer  channels 
there  was  only  one  pole  in  the  left  half  plane  aside  from  the  pole  at  the 
origin.  It  is  not  known  at  this  point  if  this  type  of  pole  configuration  will 
stay  the  same  for  arbitrary  correlation.  However,  for  small  correlation  at 
least,  one  may  expect  only  one  left  half  plane  pole.  If  this  pole  is  denoted 
by 

s  =  -I/O  (8-111) 

then ,  from  (8-110) , 

S  =  -1/a 
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Alternatively  if  the  eigenvalues  of  M  ^  are  (x^  ;  i=0,l . L)  with 

Xq=o,  then 

where  we  have  shown  explicitly  the  dependence  of  a  and  Xj^  on  x. 

It  follows  from  (8-79)  that  the  d.f.  of  p  can  be  taken  as  either  of 
the  two  equations: 


F(x) 


L 

n 

1=1 


o(x)  .-px/o(x) 
o(xy  r x'(x)  " 


F(x) 


j - sa(x) 

Jdet  U  +  sM 


g-px/o(x) 
s  =  -l/p(x) 


(8-114) 


(8-115) 


A  single  integral  expression  for  is  obtained  by  using  either 
(8-114)  or  (8-115)  in  (8-73),  with  appropriate  specification  of  (-f(p))  for 
the  modem  of  interest. 


We  will  now  examine  the  case  of  a  single  jammer  in  more  detail.  In 

this  case 


M 


01 


01 


(8-116) 


Q  = 


0 


0 


-xn. 


(8-117) 
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8.7  SOME  EXAMPLES 


For  purposes  of  illustration  we  present  some  numerical  results  for 
the  simplest  case  -  a  single  jammer  and  a  single  auxiliary  antenna.  We  exam¬ 
ine  the  case  of  PSK  transmission.  From  Eq.(4-24)  the  error  rate  is  given  by 


(8-126) 


To  evaluate  P^  it  is  necessary  to  determine  P  and  R.  For  the  case  of  one 
antenna  and  one  jammer,  these  are  given  by  (see  Eqs.  (4-4) , (4-5) , (3-4)  and 
(3-6) 


P 


1  +  IPj  1’ 


(8-127) 


(8-128) 


where  p^  is  the  optimum  complex  weight,  Sg  is  the  value  of  average  energy/bit 
/noise  power  density  for  the  received  signal,  N  is  the  thermal  noise  power  of 
the  main  and  auxiliary  antenna  outputs  and  is  the  transmittted  jamming 
power.  h^(t)  is  the  complex  gain  of  the  residual  channel  from  the  jammer  to 
the  canceller  output,  (Eq.  (2-3)) 

hi(t)  =  h^g(t)  -  Pjh^(t)  (8-129) 

in  which  h^g(t),  h^j(t)  are  the  jammer  to  main  and  auxiliary  antenna  channel 
complex  gains  (Eq.(2-6))  respectively. 
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The  optimum  weight  is  obtained  from  equations  in  Section  3.  For 


the  single  weight  case  we  find  that 


Jj  hjQh*j 


Ihjjl*  +  N 


(8-130) 


Using  (8-130)  in  (8-129)  and  computing  fh^ I '  we  find  that 


^  - 


IhiohJil 

IT  IT  ,  .  ^ 

L  IT - IT  ,  .  1 

(8-131) 


Using  the  definitions 

GlO  =  IAq(8p<'j)1* 

Gii  =  |A^(e^,«l)|* 


P  =  Sio^ll 


(8-132) 

(8-133) 

(8-134) 


for  the  antenna  gains  in  the  direction  of  the  jammer  and  using  Eq. (8-131)  and 
(8-132) 


^10^11 


|p|*  J 


10 


GjjJj  +  N 


'  ('  *  ’  ") 


(8-135) 
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Three  parameters  of  interest  are  + 


^  =  N 

(8-136) 

So 

(8-137) 

«  -  1  -  IM* 

(8-138) 

where  6  is  the  ratio  of  the  received  jammer  power  to  thermal  noise  power  in 
the  auxiliary  antenna  output,  y  is  the  ratio  of  the  main  beam  antenna  gain  to 
the  auxiliary  antenna  gain  in  the  direction  of  the  jammer,  and  6  is  a 
decorrelation  parameter. 

Using  these  parameter  definitions  in  (8-135)  and  (8-128)  we  find 

that 

R  =  ■—  (1  +  «p)  (8-139) 

Also  using  these  parameter  definitions  and  (8-132)  through  (8-133)  in  (8-130) 
and  (8-127)  we  see  that 

’’ '  ^  *  ‘w’'  ’“-‘G 

There  are  now  four  parameters  that  define  the  performance  of  the  system:  Sq, 

«,  y,  and  6. 

t 

The  decorrelation  parameter  «  is  related  to  t  in  Eq.  (6-42)  by  6  =  2e+€*  w  2e 
and  to  y  in  Eq.  (6-45)  by  6  «  2y*. 


Some  limiting  cases  are  of  interest.  Consider  first  the  case  of  no 
decorrelation,  5=0,  then 


R  =  I_ 


Ten- 


;  Y  =  0 


(8-141) 


^  t  D  *  J 


;  Y  =  0 


(8-142) 


Assume  also  a  large  jammer,  i.e.,  e  >>  1,  then 


^0 


(8-143) 


1  +  Y 


(8-144) 


1  +  D  ^0  ^  ^ 


(8-145) 


Typical  values  of  y  are  of  the  order  of  1  or  less  so  that  at  large 
SNR's  (1+P)/R  is  a  large  number.  Thus  from  Eq.(4-27) 


1  +  2y 
4s„ 


;  B  »  1,  Sn  >>  1 


(8-146) 


The  advisability  of  keeping  y  small  is  evident  from  (8-146),  the  smaller  y 
yields  the  smaller  error  rate.  It  is  also  evident  that  with  G=0  and  y  small, 
the  jammer  is  effectively  cancelled  and  the  same  performance  is  achieved  as 
with  no  jammer. 


Consider  now  the  limiting  case  6=1,  complete  decorrelation.  Then, 


R--^8  ;  e»  1,6  =  1 

0 


(B-147) 


To  the  extent  that  yB/s^  is  still  a  small  number  and  Sg  is  larger,  we  may 
still  use  the  approximation  in  Eq.(4-6), 

e  Hbg 

Now  we  see  that  the  jammer  can  increase  the  error  rate  considerably.  Under 
the  condition  6  =  i,  the  auxiliary  antenna  gain  goes  to  zero  (p^=0)  and  no 
sidelobe  cancellation  is  possible. 


Figure  8.2  presents  a  representati ve  calculation  of  error  rate  as  a 
function  of  6  for  6=1000,  y=1  with  a  family  parameter  Sg=3n  dB,  n=l,2,...9. 
Note  the  rapid  degradation  in  error  rates  with  6  at  the  high  SNR's.  Typically 
3  dB  of  degradation  sets  in  when  6  increases  from  0  to  .003  . 


I, 
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appendix 

PROOF  OF  A  MATRIX  IDENTITY 

We  shall  prove  a  matrix  identity  that  was  useful  in  the  text.  Given 
a  symmetric  matrix 


of  the  form 


^  ’  m,p  =  1,2 . Ml 


0  =  p*  =  R 

mp  pm  m-p 


(A-i; 


(A-2) 


i.e.,  Hermitian  symmetric  with  the  property  that  a  matrix  element  depends  for 
its  value  only  on  the  difference  between  its  indices,  then  the  following 
matrix  property  is  valid. 


"  det 


(A-3) 


To  prove  (A-3)  we  note  first  that  r^^^  can  be  partitioned  in  the 


form 


(A-4) 


where  we  have  normalized 


1 


(A-b) 


The  determinant  of  is  equal  to  the  sum  of  the  products  of  the 
elements  of  any  single  row,  such  as  the  first  row,  by  their  cofactors.  If  the 
cofactors  of  the  first  row  of  are  denoted  by  ••  m+1^* 

see  that 


*11  =  '•« 


(A-6) 


Pfl^l***  ^-2  ^m  ’  3  <  m  <  M-1  (A-7) 


^12  "  -215—2  ***  -^ 


(A-8) 


1,M+1 


(-1)  det  pjjj  £2  ••• 


{A-9) 


where  Cp  is  the  column  vector  in  i.e.. 


C 1  C 2  ^ 


(A-10) 


Note  that  the  matrix  in  (A-7)  differs  from  that  in  (A-8)  in  that 
the  m^^  column  of  r^  has  been  removed  and  a  new  first  column  p*  inserted.  If 

n  M 

two  columns  of  a  matrix  are  interchanged,  the  determinant  of  the  new  matrix  is 
the  negative  of  the  determinant  of  the  original  matrix.  Thus  A^^  can  be  ex¬ 
pressed  in  the  alternate  form 


|_-1’  -2‘’*  ^-2  ^  21  ’  3  £  tn  £  M  (A-11) 

^l,m+l  "  ■  j2^-l*  ^-2***  ^M-1’^  J  (A-12) 

where  we  have  interchanged  the  first  column  in  (A-7)  with  the  second  column, 

and  then  the  second  column  with  the  third,  etc.,  until  ^  interchanges  posi- 

”rl 

2 

tion  with  ^  i.e.,  m-2  interchanges  for  a  factor  of  (-1)""  . 


If  we  define  the  vector 


B 


(A-13) 


and  evaluate  the  determinant  of  r^^^  by  use  of  the  cofactors  of  the  first  row, 
it  follows  that 

det  r^^j  =  det  £^,1  (A-14) 


Examination  of  the  structure  of  I  Aj^^^  ;  m=2,...,M+l  )  reminds  one 
immediately  of  a  rule  for  the  solution  of  a  set  of  linear  equations.  In  par¬ 
ticular,  consider  the  equation  for  the  unknown  vector  x: 


Ifli  = 


(A-15) 


The  solution  for  x  is  given  by 


(A-16) 
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We  denote  as  the  p  component  of 


(A-17) 


The  solution  for  is  given  by 


det  ^+1 

det  r. 


(A-18) 


Comparison  of  (A-11)  and  (A-13)  with  (A-17)  and  (A-18)  shows  that 


^  =  X  det  ^ 


(A-19) 


where  we  have  used  the  fact  that  r^^  is  Hermitian  symmetric. 


Then,  using  (A-19)  and  (A-14)  and  dividing  the  equation  by  det  r 


we  obtain  Eq.( A-3) . 
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